A Measurement of Atmospheric Neutrino Oscillation Parameters by Super- Kamiokande I 



O 



in 



Y.Ashie/ J.Hosaka/ K.Ishihara,' Y.Itow,' J.Kameda,' Y.Koshio/ A.Minamino,' C.Mitsuda,' M.Miura,' S.Moriyama,' 
M.Nakahata,' T.Namba,' R.Nambu,^ Y.Obayashi,' M.Shiozawa,^ Y.Suzuki,' Y.Takeuchi,' K.Taki,' S.Yamada,' M.Ishitsuka,^ 
T.Kajita,2 K.Kaneyuki,^ S.Nakayama,^ A.Okada,^ K.Okumura,^ C.Saji,^ Y.Takenaga,^ S.T.Clark,^ S.Desai,^ E.Kearns,^ 
S.Likhoded,-^ J.L.Stone,^ L.R.Sulak,-^ W.Wangr^ M.Goldhaber,'^ D.Casper,^ J.P.Cravens,^ W.Gajewski,^ W.R.Kropp,^ 

D.W.Liu,^ S.Mine,^ M.B.Smy,^ H.W.Sobel,^ C.W.Stei-ner,^ M.R.Vagins,^ K.S.Ganezer,*^ J.Hill,*' W.E.Keig,*' 
J.S.JangJ J.Y.KimJ I.T.LimJ K.Scholberg,** C.W. Walter,'* R.W.Ellsworth,'' S.Tasaka,'" G.Guillian," A.Kibayashi,'' 
J.G.Learned,'' S.Matsuno,'' D.Takemori,'' M.D.Messier,'^ YHayato,'^ A.K.Ichikawa," T.Ishida,'^ T.Ishii,'-' 
T.Iwashita,'^ T.Kobayashi,'^ T.Maruyama,"'BK.Nakamura,'^ K.Nitta,'^ YOyama," M.Sakuda,"'|l| Y.Totsuka," 
A.T.Suzuki,''* M.Hasegawa,'^ K.Hayashi,'^ I.Kato,'^ H.Maesaka,'^ T.Morita,'^ T.Nakaya,'^ K.Nishikawa,'^ T.Sasaki,'^ 
S.Ueda,'^ S.Yamamoto,'^ T.J.Haines, '''•^ S.Dazeley,'^ S.Hatakeyama," R.Svoboda,'^ E.Blaufuss,'** J.A.Goodman,'** 
G.W.Sullivan,'** D.Turcan,'** A.Habig,'^ Y.Fukuda,^^ C.K.Jung,^' T.Kato,^' K.Kobayashi,^' M.Malek,^' C.Mauger,^' 
C.McGrew,2' A.SaiTat,^' E.Sharkey,^' C.Yanagisawa,^' T.Toshito,^^ K.Miyano,^^ N.Tamura,^^ J.Ishii,24 Y.Kuno,24 
O ■ M.Yoshida,24 S.B.Kim,^^ J.Yoo,^^ H.Okazawa,^^ T.Ishizuka,^^ Y.Choi,^^ H.K.Seo,^** Y.Gando,^'^ T.Hasegawa,^'' 

CN ! K.Inoue,29 J.Shirai,29 A.Suzuki,^'' M.Koshiba,^'' Y.Nakajima,^' K.Nishijima,^' T.Harada,^^ H.Ishino,^^ Y.Watanabe,^^ 

pH ■ D.Kielczewska,^^'^ J.Zalipska,^^ H.G.Berns,34 R.Gran,34 K.K.Shiraishi,34 A.Stachyra,^^ K.Washburn,34 and R.J.Wilkes^^ 
^ ^ ' (The Super-Kamiokande Collaboration) 

^ Kamioka Observatory, Institute for Cosmic Ray Research, University of Tokyo, Kamioka, Gifu, 506-1205, Japan 
^ ^Research Center for Cosmic Neutrinos, Institute for Cosmic Ray Research, University of Tokyo, Kashiwa, Chiba 277-8582, Japan 

^Department of Physics, Boston University, Boston, MA 02215, USA 
' '^Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA 

^ ' ^ Department of Physics and Astronomy, University of California, Irvine, Irvine, CA 92697-4575, USA 

, ^Department of Physics, California State University, Dominguez Hills, Carson, CA 90747, USA 

' ^Department of Physics, Chonnam National University, Kwangju 500-757, Korea 

I ^Department of Physics, Duke University, Durham, NC 27708, USA 

. ^Department of Physics, George Mason University, Fairfax, VA 22030, USA 

' '^Department of Physics, Gifu University, Gifu, Gifu 501-1193, Japan 

, Department of Physics and Astronomy, University of Hawaii, Honolulu, HI 96822, USA 

' '^Department of Physics, Indiana University, Bloomington, IN 47405-7105, USA 

^ ' High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan 

Department of Physics, Kobe University, Kobe, Hyogo 657-8501, Japan 
^ ' '^Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

^ I '^Physics Division, P-23, Los Alamos National Laboratory, Los Alamos, NM 87544, USA 

^ ■ '^Department of Physics and Astronomy, Louisiana State University, Baton Rouge, LA 70803, USA 

^ ' '^Department of Physics, University of Maryland, College Park, MD 20742, USA 

• , '^Department of Physics, University of Minnesota, Duluth, MN 55812-2496, USA 

' ^'^ Department of Physics, Miyagi University of Education, Sendai,Miyagi 980-0845, Japan 

5— ( ' ^'Department of Physics and Astronomy, State University of New York, Stony Brook, NY 11794-3800, USA 

5^ , Department of Physics, Nagoya University, Nagoya, Aichi 464-8602, Japan 

Department of Physics, Niigata University, Niigata, Niigata 950-2181, Japan 
Department of Physics, Osaka University, Toyonaka, Osaka 560-0043, Japan 
Department of Physics, Seoul National University, Seoul 151-742, Korea 
^^International and Cultural Studies, Shizuoka Seika College, Yaizu, Shizuoka 425-8611, Japan 
Department of Systems Engineering, Shizuoka University, Hamamatsu, Shizuoka 432-8561, Japan 
Department of Physics, Sungkyunkwan University, Suwon 440-746, Korea 
^'^ Research Center for Neutrino Science, Tohoku University, Sendai, Miyagi 980-8578, Japan 
^'^ University of Tokyo, Tokyo 113-0033, Japan 
^'Department of Physics, Tokai University, Hiratsuka, Kanagawa 259-1292, Japan 
Department of Physics, Tokyo Institute for Technology, Meguro, Tokyo 152-8551, Japan 
Institute of Experimental Physics, Warsaw University, 00-681 Warsaw, Poland 
^''Department of Physics, University of Washington, Seattle, WA 98195-1560, USA 
(Dated: February 3, 2008) 

We present a combined analysis of fully-contained, partially-contained and upward-going muon atmospheric 
neutrino data from a 1489 day exposure of the Super-Kamiokande detector. The data samples span roughly 
five decades in neutrino energy, from 100 MeV to 10 TeV. A detailed Monte Carlo comparison is described and 
presented. The data is fit to the Monte Carlo expectation, and is found to be consistent with neutrino oscillations 
of ^ Vt with sin22e > 0.92 and 1.5 x 10"^ < Am^ < 3.4 x lO^^eV^ at 90% confidence level. 
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PACS numbers: PACS numbers: 14.60.Pq, 96.40.Tv 



I. INTRODUCTION 

Atmospheric neutrinos are produced from the decays of 
particles resulting from interactions of cosmic rays with 
Earth's atmosphere. We have previously reported the results 
of a number of atmospheric neutrino observations spanning 
energies from 100 MeV to 10 TeV yj,i2|,i3|,j4J. In each case, a 
significant zenith-angle dependent deficit of was observed. 
These deficits have been interpreted as evidence for neutrinos 
oscillations |5]. If neutrinos have a non-zero mass, then the 
probability that a neutrino of energy produced in a weak 
flavor eigenstate Va will be observed in eigenstate Vp after 
traveling a distance L through the vacuum is: 



^'(Va ^Vp) = sin2 2esin^ 



1.27Am2(eV2)L(km)\ 



£v(GeV) 



/ 7 



(1) 



where 9 is the mixing angle between the mass eigenstates 
and the weak eigenstates and Am^ is the difference of the 
squared mass eigenvalues. This equation is valid in the 2- 
flavor approximation. The analysis reported in this paper is 
under the assumption of effective 2-flavor neutrino oscilla- 
tions, <-> Vx, which is considered to be dominant in atmo- 
spheric neutrino oscillations. Equation^is also true in matter 
for Vx, but may be modified for oscillation involving Vg 
which travel through matter The zenith angle dependence of 
the observed deficits results from the variation of L with the 
direction of the neutrino. Neutrinos produced directly over- 
head travel roughly 15 km to the detector while those pro- 
duced directly below traverse the full diameter of the Earth 
(13,000 km) before reaching the detector By measuring the 
neutrino event rate over these wide ranges of and L, we 
have measured the neutrino oscillation parameters Am^ and 

Super-Kamiokande (also Super-K or SK) is a 50-kiloton 
water Cherenkov detector located deep underground in Gifu 
Prefecture, Japan. Atmospheric neutrinos are observed in 
Super-K in two ways. At the lowest energies, 100 MeV - 
10 GeV, atmospheric neutrinos are observed via their charged- 
current interactions with nuclei in the 22.5 kiloton water fidu- 
cial mass: v + N ^ I + X. These interactions are classified 
as fully-contained (FC) if all of the energy is deposited inside 
the inner Super-K detector, or as partially-contained (PC) if a 
high energy muon exits the inner detector, depositing energy 
in the outer veto region. The neutrino energies that produce 
partially-contained events are typically 10 times higher than 
those that produce fully-contained events. The Super-K detec- 
tor started observation on April, 1996 achieving a 92kiloton- 
yr (1489 live-day) exposure to atmospheric neutrinos through 
July, 2001 during the Super-Kamiokande 1 running period. 
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FIG. 1: The parent neutrino energy distributions for the fully- 
contained, partially-contained, upward stopping-muon and upward 
through-going muons samples. Rates for the fully-contained and 
partially-contained samples are for interactions in the 22.5 kiloton 
fiducial volume. Taken together, the samples span five decades in 
neutrino energy. 



Neutrinos can also be detected by their interactions with 
the rock surrounding the detector. Charged-current in- 
teractions with the rock produce high energy muons which 
intersect the detector While these interactions can not be 
distinguished from the constant rain of cosmic ray muons 
traveling in the downward direction, muons traveling in an 
upward direction through the detector must be neutrino in- 
duced. Upward-going muon events are separated into two cat- 
egories: those that come to rest in the detector (upward stop- 
ping muons) and those that traverse the entire detector volume 
(upward through-going muons). The energies of the neutrinos 
which produce stopping muons are roughly the same as for 
partially-contained events, ~ 10 GeV. Upward through-going 
events, however, are significantly more energetic; the parent 
neutrino energy for these events is about 100 GeV on average. 

Figure n shows the expected number of neutrino events in 
each event category as a function of neutrino energy. The 
samples taken together span nearly five decades in energy. 
This broad range of available energies, in combination with 
the variation in neutrino travel distance, makes the combined 
data sample well suited for a precise measurement of neutrino 
oscillation parameters. 

There have been numerous other measurements of atmo- 
spheric neutrinos. Kamiokande 't'], 1MB fsl ^ and 
Soudan 2 1 10, 11] observed significantly smaller to Vg flux 
ratios of ^ 1 GeV atmospheric neutrinos, which were inter- 
preted as a signature for neutrino oscillation. The ratio was 
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used in order to normalize the uncertainty in the overall atmo- 
spheric neutrino flux. Data on multi-GeV atmospheric neu- 
trino events fl2ll and upward-going muons llisl Il4l ITsIl have 
also shown a zenith-angle dependent deficit of the flux. The 

^ Vx oscillation analyses of these various data over various 
energy ranges OHBIlSllllilEHIil indicated sim- 
ilar Am^ and sin^ 29 regions as the first measurements from 
Super-K as well as those reported here. 

The K2K long baseline experiment used an accelerator 
beamline to produce muon neutrinos that traveled 250 km to 
the Super-K detector, as a means to study neutrino oscillation 
in the atmospheric neutrino energy and distance scales. The 
results from K2K 1 19, 20] are also consistent with the neutrino 
oscillation parameters reported here. 



II. ATMOSPHERIC NEUTRINOS 

To carry out detailed studies of neutrino oscillations using 
atmospheric neutrinos, it is important to know the expected 
flux without neutrino oscillations. The difficulties and the un- 
certainties in the calculation of atmospheric neutrino fluxes 
differ between high and low energies. For low energy neutri- 
nos around 1 GeV, the primary fluxes of cosmic ray compo- 
nents are relatively well known. Low energy cosmic ray fluxes 
of less than about 10 GeV are modulated by solar activity, with 
the minimum flux occurring at times of high solar activity. At 
these energies, the primary cosmic rays are also affected by 
the geomagnetic field through a rigidity (momentum/charge) 
cutoff. For high energy neutrinos, above 100 GeV, primary 
cosmic rays with energies greater than 1000 GeV are relevant. 
At these energies, solar activity and the rigidity cutoff do not 
affect the cosmic rays, but details of the higher energy primary 
cosmic ray flux are not as well measured. 

There are several flux calculations |21, 22, 23, 24, 25, '26', 
mmiH. Unlike older calculations |3CL,HJ, in which the 
secondary particles were assumed to travel in the direction 
of the primary cosmic ray (1 -dimensional calculations), the 
current calculations employ three dimensional Monte Carlo 
methods. We outline below the methods of the calculation. 
We compared results from three atmospheric neutrino flux 
calculations |25, 28, 29] which cover the energy range rele- 
vant to the present analysis. The flux from Honda et al. ] 28] is 
used for the main numbers and figures quoted for the Super- 
Kamiokande analysis. 

Calculations start with primary cosmic rays based on mea- 
sured fluxes, and include solar modulation and geomagnetic 
field effects. The interaction of cosmic ray particles with the 
air nucleus, the propagation and decay of secondary particles 
are simulated. We used a neutrino flux calculated specifi- 
cally for the Kamioka site. According to the cosmic ray pro- 
ton, helium and neutron measurements the cosmic 
ray flux was near that of solar minimum until the summer of 
1999, rapidly decreased during the next year, and was at the 
minimum value consistent with solar maximum from summer 
of 2000 until Super-Kamiokande stopped taking date in July 
2001 . Therefore, the atmospheric neutrino Monte Carlo is cal- 
culated for 3 years of solar minimum, 1 year of changing ac- 



tivity, and 1 year of solar maximum. 
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FIG. 2: (a) The direction averaged atmospheric neutrino energy spec- 
trum for + V^ calculated by several authors are shown by solid 
line 1 28], dashed line |29] and dotted line 13] ■ © The ratio of the 
calculated neutrino flux. The fluxes calculated in r29] (solid line) and 
|23l (dashed line) are normalized by the flux in |28]. 

The calculated energy spectra of atmospheric neutrinos at 
Kamioka are shown in Fig. |2ja). Also shown in Fig. |2jb) is 
the comparison of the calculated fluxes as a function of neu- 
trino energy. The agreement among the calculations is about 
10 % below 10 GeV. This can be understood because the accu- 
racy in recent primary cosmic ray flux measurements ] 34, 3^ 
below 100 GeV is about 5 % and because hadronic interaction 
models used in each calculation are different. 

However the primary cosmic ray data are much less accu- 
rate above 100 GeV. Therefore, for neutrino energies much 
higher than 10 GeV, the uncertainties in the absolute neutrino 
flux could be substantially larger than the disagreement level 
among the calculations. In Ref. ]36], the authors discussed 
that the fit to the low energy (< 100 GeV) proton spectra gave 
a spectrum index of -2. 74-1-0.01 . However, this spectrum does 
not fit well to the high energy data. Therefore, authors in 
Ref. ]|28] fit the high energy data allowing a different spectral 
index above 100 GeV and found the best fit value of -2.71. 
There is 0.03 difference in the spectrum index for low energy 
(<100 GeV) and high energy (> 100 GeV) protons. Also, it is 
discussed in Ref. ]36] that the spectrum index for the He flux 
can be fit by either -2.64 or -2.74. There could be 0.10 uncer- 
tainty in the spectrum index for He. The spectrum indices for 
heavier nuclei have uncertainties larger than 0.05 ]36]. Taking 
the flux weighted average of these spectrum index uncertain- 
ties, we assign 0.03 and 0.05 for the uncertainties in the energy 
spectrum index in the primary cosmic ray energy spectrum be- 
low and above 100 GeV, respectively. 

Figure|3lshows the calculated flux ratio of + to Ve + Vg 
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as a function of the neutrino energy, integrated over solid an- 
gle. This ratio is essentially independent of the primary cos- 
mic ray spectrum. Especially in the neutrino energy region of 
less than about 5 GeV, most of the neutrinos are produced by 
the decay chain of pions and the uncertainty of this ratio is 
about 3 %, which is estimated by comparing the three calcula- 
tion results. The contribution of K decay in neutrino produc- 
tion is more important in the higher energy region; about 10 % 
for + Vg and 20 % for v^, + Vu at 1 GeV. It increases to more 
than 30 % at 100 GeV for both + Vg and + V^,. There, the 
ratio depends more on the K production cross sections and the 
uncertainty of the ratio is expected to be larger. A 20 % un- 
certainty in the K/ti production ratio I3ll 13711 causes at least 
a few percent uncertainty in the + to Vg + Vg ratio in the 
energy range of 10 to 100 GeV. However, as seen from Fig.|3j 
the difference in the calculated v^, + Vu to Vg + ratio is as 
large as 10% at 100 GeV. As a consequence, above 5 GeV, 
we assumed that the uncertainty linearly increases with logEv 
from 3 % at 5 GeV to 10 % at 100 GeV. 
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FIG. 3: The flux ratio of + to Vg + Vg averaged over all zenith 
and azimuth angles versus neutrino energy. Solid, dashed and dotted 
lines show the prediction by Ii28<1 . L29.1 and 1.25.1 . respectively. 



Figure |4] shows the calculated flux ratios of to and Vg 
to Vg. The calculations agree to about 5 % for both of these 
ratios below 10 GeV. However, the disagreement gets larger 
above 10 GeV as a function of neutrino energy. The sys- 
tematic errors in the v/v ratio are assumed to be 5 % below 
10 GeV and linearly increase with log/iy to 10 % and 25 % at 
100 GeV, for the Vg to Vg and v^j to ratios, respectively. 

Figure |5] shows the zenith angle dependence of the at- 
mospheric neutrino fluxes for several neutrino energies. At 
low energies, and at the Kamioka location, the fluxes of 
downward-going neutrinos are lower than those of upward- 
going neutrinos. This is due to the deflection of primary cos- 
mic rays by the geomagnetic field, roughly characterized by a 
minimum rigidity cutoff. For neutrino energies higher than a 
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FIG. 4: The flux ratios of to and Vg to Vg versus neutrino energy. 
Solid, dashed and dotted lines show the prediction by |2^ . |2^ and 
l23l . respectively (same key as Fig.|3j. 



few GeV, the calculated fluxes are essentially up-down sym- 
metric, because the primary particles are more energetic than 
the rigidity cutoff. 

The enhancement of the flux near horizon for low energy 
neutrinos is a feature characteristic of the three dimensional 
nature of the neutrino production in cosmic ray hadronic 
showers. This is properly treated in current flux calcula- 
tions liT, 22, 23, 24, 25, "2?, Uf, "28", "291. However, in 
Super-Kamiokande, the horizontal enhancement cannot be 
seen in the lepton zenith angle distribution, due to the rela- 
tively poor angular coiTelation between neutrinos and leptons 
below 1 GeV. The uncertainties in the up-down and vertical- 
horizontal ratios of the number of events are estimated by 
comparing the predicted ratios by various flux models. These 
uncertainties generally depend on the energy and the neutrino 
flavor. The uncertainty in the up-down event ratio is about 1 
to 2 % in the sub-GeV energy region and is about 1 % in the 
multi-GeV energy region. The main source of the uncertainty 
in the vertical-horizontal ratio around a GeV is the size of the 
horizontal enhancement of the flux due to the three dimen- 
sional effect; the uncertainty is estimated to be less than a few 
percent. 

In the higher energy region, where upward through-going 
muons are relevant, the largest source of the uncertainty in 
the vertical-horizontal ratio is the K production cross section. 
We assume that the K/n production ratio uncertainty is 20 % 
in the whole energy region IJH l3^ . The uncertainties in the 
zenith angle and energy distributions due to the K/n produc- 
tion uncertainty are included in the systematic errors in the 
analysis. This error is most important for higher energy neu- 
trinos. For example, the vertical-horizontal uncertainty for 
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upward through-going muons due to the K/ti production un- 
certainty is estimated to be 3 % L38.1 . Figure |5] shows the 
zenith angle dependence of the atmospheric neutrino fluxes 
for higher energy region observed as upward muons in Super- 
Kami okande. 
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FIG. 5: The flux of atmospheric neutrinos versus zenith angle. Solid, 
dashed and dotted lines show the prediction by 12^1 . and |2^ . 
respectively (same key as Fig.|3j- 




FIG. 6: The flux of upward-going atmospheric neutrinos versus 
zenith angle for higher energy region. Solid, dashed and dotted lines 
show the prediction by i28il . l29ll and l2^ . respectively (same key as 
Fig.|3j. 



The flight length of neutrinos is an important ingredient in 
the analysis of neutrino oscillation. For neutrinos passing a 
great distance through the Earth, the flight length can be ac- 
curately estimated. However, for horizontal and downward 
going neutrinos, the height of production in the upper atmo- 
sphere must be distributed by the Monte Carlo method. Fig- 
ure0shows the calculated flight length distributions for verti- 
cally down-going and horizontal neutrinos. 

In summary of the atmospheric neutrino flux, we remark 
that, the (v^^, + Vu) over (Ve + V,.) flux ratio is predicted to an 
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FIG. 7: The calculated flight length of neutrinos for vertically down- 
going (cos© = 0.95 — 1.00) and near horizontal-going (cos© = 
0.05 — 0.10) directions. Distributions are made for both muon- 
neutrinos and electron-neutrinos and for three energy intervals [28]. 



accuracy of about 3 % in the energy region relevant to the data 
analysis discussed in this paper The zenith angle dependence 
of the flux is well understood, and especially, above a few 
GeV neutrino energies, the flux is predicted to be up-down 
symmetric. 



III. THE SUPER-KAMIOKANDE DETECTOR 

Super-Kamiokande is a 50 kiloton water Cherenkov detec- 
tor located at the Kamioka Observatory of the Institute for 
Cosmic Ray Research, University of Tokyo. Figure|8]shows a 
cut-away diagram of the Super-Kamiokande detector This fa- 
cility is in the Mozumi mine of the Kamioka Mining Company 
in Gifu prefecture, Japan, under the peak of Mt. Ikenoyama, 
providing a rock overburden of 2,700 m.w.e. Super-K con- 
sists of two concentric, optically separated water Cherenkov 
detectors contained in a stainless steel tank 42 meters high and 
39.3 meters in diameter, holding a total mass of 50,000 tons of 
water. The inner detector (ID) is comprised of 11, 146 Hama- 
matsu R3600 50 cm diameter photomultiplier tubes (PMTs), 
viewing a cylindrical volume of pure water 16.9 m in ra- 
dius and 36.2 m high. The 50 cm PMTs were specially de- 
signed f3^ to have good single photoelectron (p.e.) response, 
with a timing resolution of 2.5 nsec RMS. The ID is sur- 
rounded by the outer detector (OD), a cylindrical shell of wa- 
ter 2.6 to 2.75 m thick including a dead space 55 cm. The 
OD is optically isolated from the ID, and is instrumented with 
1,885 outward-facing Hamamatsu R1408 20 cm PMTs, pro- 
viding both a 47t active veto and a thick passive radioactiv- 
ity shield. The information from the outer detector is used to 
identify both incoming and outgoing muons. 

Both ID and OD PMT signals are processed by asyn- 



6 




FIG. 8: A drawing of the Super- Kamiokande detector The cutaway 
shows the inside hned with photomultiplier tubes comprising a pho- 
tocathode coverage of about 40%. The support structure is stainless 
steel beams. The thin outer region is shown with sparser density of 
outward facing PMTs: 2 outer PMTs for every unit of 3 x 4 inner 
PMTs. The top of the detector, under the hemisphere, consists of 
electronics huts and open work area. 



chronous, self-triggering circuits that record the time and 
charge of each PMT hit over a threshold. Each ID PMT signal 
is digitized with custom Analog Timing Modules (ATMs) |40, 
14111 which provide 1 .2 /jsec timing range at 0.3 nsec resolution 
and 550 pC charge range at 0.2 pC resolution 0.1 p.e.). 
The ATM has automatically-switched dual channels to pro- 
vide deadtime-free data acquisition. The outer PMT signals 
are processed with custom charge-to-time conversion mod- 
ules and digitized with LeCroy 1877 multi-hit TDCs over a 
— 10 /jsec to +6 /jsec window centered on the trigger time. 
More details of the Super-K detector can be found in |42|. 

An event used in the atmospheric neutrino analysis is trig- 
gered by the coincidence of at least 30 PMT hits in a 200 nsec 
window. The hit threshold for each individual PMT is about 
1 /4 p.e. This trigger condition corresponds to the mean num- 
ber of hit PMTs for a 5.7 MeV electron. The trigger rate is 
10-12 Hz. The trigger rate due to cosmic ray muons is 2.2 Hz. 
Digitized data are saved at a total rate of 12 GB per day. 

The detector is simulated with a Monte Carlo program 
based on the GEANT package i43ll . in which the propaga- 
tion of particles, the generation and propagation of Cherenkov 
photons, and the response of the PMTs is considered. For 
hadronic interactions in water, the CALOR package ll44ll was 
employed in our simulation code. This package is known 
to reproduce the pion interactions well down to low mo- 
mentum regions of 1 GeV/c. For still lower momenta 
iPn < 500 Me V/c), a custom program based on experimental 



data from 7t — '^O scattering i45ll and % — p scattering ll46ll was 
used in our simulation code. 

In connection with the propagation of charged particles, 
Cherenkov photons are generated. For the propagation of 
Cherenkov photons in water, Rayleigh scattering, Mie scat- 
tering and absorption were considered in our simulation code. 
The attenuation coefficients used were tuned to reproduce the 
measurement using laser system (Section flll A> . Light reflec- 
tion and absorption on detector material, such as the surface 
of PMTs and black plastic sheets between the PMTs was sim- 
ulated based on direct measurements, using probability func- 
tions that depend on the photon incident angle. 



A. Calibration of the Super-Kamiokande Detector 

Water transparency was measured using a dye laser beam 
injected into detector water at wavelengths of 337, 371, 400, 
and 420 nm. From the spatial and timing distribution of ob- 
served laser light, both absorption and scattering coefficients 
were studied and incorporated into our detector simulator. The 
water transparency was continuously monitored using cosmic 
ray muons as a calibration source. 

The accuracy of the absolute energy scale was estimated 
to be ±1.8% based on the following calibration sources: the 
total number of photo-electrons as a function of muon track 
length, where the muon track length is estimated by the re- 
constructed muon entrance point and the reconstructed vertex 
point of an electron from the muon-decay; the total number 
of photo-electrons as a function of Cherenkov angle for low 
energy cosmic ray muons; the spectrum of muon-decay elec- 
trons; and the invariant mass of T\Ps produced by neutrino in- 
teractions (Figure|9}. Figure [Tol summarizes the absolute en- 
ergy scale calibration by these studies. The stability of the 
energy scale was also monitored continuously using stopping 
muons and muon-decay electrons. Figure ^2 shows the time 
variation of the mean reconstructed energy of stopping muons 
divided by muon range and the mean reconstructed electron 
energy from muon-decays. The R.M.S of the energy scale 
variation is ±0.9% over the time of the experiment. From 
combining the absolute energy scale accuracy study (±1 .8 %) 
and the energy scale time variation (±0.9%), the total uncer- 
tainty of the energy scale of atmospheric neutrino detection 
was estimated to be ±2.0 %. 

The uniformity of the detector response was studied by 
decay electrons from stopping cosmic ray muons and neu- 
trino induced Tt". Both are good calibration sources because 
the vertex position is distributed in the fiducial volume and 
the momentum distribution is nearly uniform in all direc- 
tions. To account for muon polarization in the estimation 
of the zenith and azimuthal angle dependence of the detec- 
tor gain, only electrons decaying in the direction perpendic- 
ular to the initial muon direction are used. This condition is 
—0.25 < cos0j,M^ < 0.25 where ©g^^ is the opening angle 
between the electron and muon directions. Using the selected 
electrons, the mean of the reconstructed momentum of the 
electrons are plotted as a function of the zenith angle of the 
electrons in Figure [T2l-(a). From the fi gure, the detector gain 
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FIG. 11: The mean reconstructed energy of cosmic ray stopping 
muons divided by their range (upper) and muon-decay electron 
(lower) as a function of elapsed days. Vertical axes in both figures 
are normalized to mean values and each data point corresponds to 
two month period. The variation is within ± 2 %. 
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FIG. 9: The invariant mass distribution of fully-contained events with 
two e-like rings and no muon-decay electron, for SK data (points) 
and atmospheric neutrino Monte Carlo (histogram). A peak from 
neutrino induced 7l" is clearly observed. 
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FIG. 10: The determination of the absolute energy scale of Super- 
Kamiokande based on in situ calibration with jU-decay electrons, 
TjO _> yy invariant mass, and the Cherenkov light of stopping cosmic 
ray muons. 
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FIG. 12: The gain uniformity of the Super-Kamiokande detector as 
determined by the mean value of the reconstructed decay electron 
momentum (a) as a function of zenith angle, and (b) as a function of 
azimuthal angle. Vertical axes in both figures are normalized to the 
mean values. 



IV. ATMOSPHERIC NEUTRINO MONTE CARLO 



was uniform over all zenith angles within ±0.6 %. Figure 
I12l-(b) shows the azimuthal angle dependence of the recon- 
structed momentum. Again, the detector gain is uniform over 
all azimuthal angles within ±1 %. Finally, Figure [O] shows 
the zenith angle dependence of the reconstructed 71'' mass. 
This figure also suggests that the detector gain was uniform 
over all zenith angles within ±1 %. 



The result published in this paper relies heavily on detailed 
comparison of the experimental data with the theoretical ex- 
pectation. An important element of this is to simulate the 
interaction of neutrinos from 10 MeV to 100 TeV with the 
nuclei of water, or in the case of upward muons, the nuclei 
of the rock surrounding the detector, assumed to be "stan- 
dard rock"(Z=ll, A=22). Therefore, we have developed two 
Monte Carlo models designed to simulate neutrino interac- 
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FIG. 13: The gain uniformity of the Super- Kamiokande detector as 
determined by the fitted peak of the Jt" ^ YY mass distribution as 
a function of zenith angle. Vertical axis is normalized to the mean 
value. 



tions with protons, oxygen and sodium 147114811 . Both models 
use similar input physics models. Here, one of the models il 
(NEUT) will be described. The detailed description for the 
other model (NUANCE) can be found elsewhere |48J. 

In the simulation program, the following charged and neu- 
tral current neutrino interactions are considered: 

• (quasi-)elastic scattering, v N ^ I N', 

• single meson production, v N ^ I N' m, 

• coherent 7i production, v '^O I K '^O, 

• deep inelastic scattering, vA^ I N' hadrons. 

Here, and A^' are the nucleons (proton or neutron), / is the 
lepton, and m is the meson, respectively. For single meson 
production, K and r| production are simulated as well as the 
dominant n production processes. If the neutrino interaction 
occurred in the oxygen nuclei, generated particles like pions 
and kaons interact with the nucleus before escaping. 

A. Elastic and quasi-elastic scattering 

The formalization of quasi-elastic scattering off a free pro- 
ton, which was used in the simulation programs, was de- 
scribed by Llewellyn-Smith]^^. For scattering off nucleons 
in '^O, the Fermi motion of the nucleons and Pauli Exclusion 
Principle were taken into account. The nucleons are treated as 
quasi-free particles using the relativistic Fermi gas model of 
Smith and MonizfsO]. The momentum distribution of the nu- 
cleons were assumed to be flat up to the fixed Fermi surface 
momentum of 225 MeV/c. This Fermi momentum distribu- 
tion was also used for other nuclear interactions. The nuclear 
potential was set to 27 MeV/c. 

B. Single meson production 

Rein and Sehgal's model was used to simulate the reso- 
nance productions of single 7t, K and r] llsil l52il . In this 
method, the interaction is separated into two parts: 



N* 



l+N*, 
m+N\ 



where m is a meson, A^ and A^' are nucleons, and A^* is a baryon 
resonance. The hadronic invariant mass, W, the mass of the 
intermediate baryon resonance, is restricted to be less than 
2 GeV/c^. In addition to the dominant single K production, 
K and r| production is considered. The production of r| is evi- 
dently much smaller than n, as seen in Fig.|5]where there is no 
evidence for a mass peak near 549 MeV/c^ in data or Monte 
Carlo. 

To determine the angular distribution of pions in the final 
state, we also use Rein and Sehgal's method for the ^33(1232) 
resonance. For the other resonances, the directional distribu- 
tion of the generated pions is set to be isotropic in the res- 
onance rest frame. The angular distribution of 7t+ has been 
measured for vp /j^pK^ \ 53] and the results agree well with 
the Monte Carlo prediction. We also consider the Pauli block- 
ing effect in the decay of the baryon resonance by requiring 
that the momentum of nucleon should be larger than the Fermi 
surface momentum. Pion-less delta decay is also considered, 
where 20 % of the events do not have the pion and only the 
lepton and nucleon are generated 13 ■ 

The quasi-elastic and single meson production models have 
a parameter (axial vector mass. Ma) that must be determined 
by experiments. For larger values, interactions with higher 
values (and therefore larger scattering angles) are en- 
hanced for these channels. The Ma value was tuned using 
the K2K fl^ near detector data. In our atmospheric neutrino 
Monte Carlo simulation. Ma is set to 1 . 1 GeV for both the 
quasi-elastic and single-meson production channels, but the 
uncertainty of the value is estimated to be 10%. Figure fT^ 
shows the K2K 1 kton water Cherenkov data on the scattering 
angle for single Cherenkov ring events fl9il together with the 
prediction by the Monte Carlo used in this analysis. The scat- 
tering angle agrees well between the data and Monte Carlo 
overall, although the suppression of events at small angle is 
being studied by several groups]55]. 

Coherent single-pion production, the interaction between 
the neutrino and the entire oxygen nucleus, is simulated using 
the formalism developed by Rein and Sehgal Il5611 . 



C. Deep inelastic scattering 

In order to calculate the cross-sections of deep inelastic 
scattering, the GRV94 IstIi parton distribution function is 
used. In the calculation, the hadronic invariant mass, W, is 
required to be greater than 1.3 GeV/c^. However, the multi- 
plicity of pions is restricted to be greater than or equal to 2 
for 1.3 < W < 2.0 GeV/c2, because single pion production 
is separately simulated as previously described. In order to 
generate events with multi-hadron final states, two models are 
used. For W between 1.3 and 2.0 GeV/c^, a custom-made 
program is used to generate the final state hadrons; only 
pions are considered in this case. For W larger than 2 GeV/c^, 
PYTHIA/JETSET ]59] is used. 

Total charged current cross sections including quasi-elastic 
scattering, single meson productions and deep inelastic scat- 
tering are shown in Figll5l 
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FIG. 14: The scattering angle distribution by neutrino interactions 
off the H^O target from the K2K experiment (data are from Fig. 1(b) 
of Ref. |T3])- Single Cherenkov ring events observed by the 1 kton 
water Cherenkov detector are used. The histogram shows the predic- 
tion by the Monte Carlo used in the present analysis. 



D. Nuclear effects 



The interactions of mesons within the '^O nucleus are also 
important for the atmospheric neutrino analysis. Basically, 
all of the interactions are treated by using a cascade model. 
The interactions of pions are very important because the cross 
section for pion production is quite large for neutrino energies 
above 1 GeV and the interaction cross sections for pions in 
nuclear matter is also large. 

In our simulation program, we consider the following pion 
interactions in '''O: inelastic scattering, charge exchange and 
absorption. The procedure to simulate these interactions is as 
follows. The initial position of the pion generated according 
to the Woods-Saxon nucleon density distribution! 72]. The in- 
teraction mode is determined from the calculated mean free 
path of each interaction. To calculate the mean free path, we 
adopt the model described by Salcedo et al. |73]. The calcu- 
lated mean free path depends not only on the momentum of 
the pion but also on the position of the pion in the nucleus. 
If inelastic scattering or charge exchange occurs, the direction 
and momentum of the pion are determined by using the re- 
sults of a phase shift analysis obtained from K — N scattering 
experiments 1 74 1 . When calculating the pion scattering ampli- 
tude, the Pauli blocking effect is also taken into account by 
requiring the nucleon momentum after interaction to be larger 
than the Fermi surface momentum at the interaction point. 
The pion interaction simulation was checked using data for 
the following three interactions: Tt^^C scattering, Tt'^O scat- 
tering and pion photo-production (y+'^C — > +X) L75il . 
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FIG. 15: Charged current total cross section divided by Ey for (a) 
neutrino and (b) anti-neutrino nucleon charged current interactions. 
Solid line shows the calculated total cross section. The dashed, dot 
and dash-dotted lines show the calculated quasi-elastic, single-meson 
and deep-inelastic scatterings, respectively. Data points are taken 
from the following experiments: (A)ANL|60], (0)GGM77f6ll], 
(•)GGM79(a)|62],(b)|63], (*)Serpukhov| 64j, (0)ANL82|6a], 
(★)BNL86|53], (■)CCFR90|66], (T)CDHSW87|6i, 

(x)IHEP-JINR96l62J, (+)IHEP-ITEP79|69|l, (□)CCFRR84|^, 
and (A)BNL82|2ir 



V. SUPER-KAMIOKANDE DATA 

The Super-Kamiokande data set was acquired from May, 
1996 to July, 2001. Three separate data reduction paths were 
used to isolate samples of fully-contained events, partially- 
contained events, and upward-going muons. The fully- 
contained and partially-contained data sets shared a common 
set of good run selection criteria, and have identical live-time. 
The upward-going muon data set relies mostly on fitting long 
muon track directions; it was less susceptible to detector ef- 
fects, and therefore had looser data quality cuts and somewhat 
higher live-time. 
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To separate fully-contained and partially-contained events, 
a fast spatial clustering algorithm was applied to the outer 
detector hits; if the number of hits in the largest OD clus- 
ter was less than 10, the event was defined as fully-contained 
(FC), otherwise, it was defined as partially-contained (PC). 
Figure ^] shows the number of the outer detector hits in the 
largest OD cluster. A clear separation of FC and PC events 
is seen at 10 hits. The systematic uncertainty of the FC and 
PC separation was estimated by scaling the number of outer 
detector hits to match the distribution among data and MC. 
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FIG. 16: The number of hits in the largest outer detector cluster, 
which is used to separate the fully-contained and partially-contained 
event samples. The histogram shows the MC prediction with neu- 
trino oscillations. 

In the early stage of the Super-Kamiokande experiment, 
two data analyses based on independent data reduction, re- 
construction, and simulation were carried out to make sure 
that the atmospheric neutrino results did not have any serious 
mistakes |1, 2]. After confirming that the two analyses de- 
scribed atmospheric neutrino results equally well, they were 
unified. This paper is based on the unified analysis plus meth- 
ods developed since that time. 



A. Fully-Contained Events 

1. Data Reduction 



Cherenkov light (for example, a muon of momentum 300 
MeV/c yields a visible energy of about 110 MeV). To make 
the final FC data sample, five steps of data reduction criteria 
were used: 

(i-ii) Simple and efficient criteria were applied in the first 
and second reduction steps: (1) the total charge collected in 
the inner detector within a 300 nsec time window must be 
greater than 200 p.e.'s; (2) the ratio of the maximum p.e. in 
any single ID PMT to the total number of p.e.'s in the inner 
detector must be less than 0.5; (3) the number of hits in the 
outer detector within an 800 nsec time window should be less 
than 25 hits; (4) the time interval from the preceding event 
should be greater than 100 /jsec, to reject electrons from stop- 
ping muon decays. 

(iii) More complex criteria were applied in the third reduc- 
tion step with the help of event reconstruction tools, for further 
rejection of cosmic ray muons and low energy events: (1) no 
spatial cluster of more than 10 OD PMT hits is allowed within 
8 m from the entrance or exit point of a candidate muon track 
fit to the inner detector light pattern, (2) the number of ID hits 
in 50 nsec residual time window should be 50 hits or more. 

(iv) In the fourth reduction step, additional selection cri- 
teria were used to eliminate spurious events, such as those 
due to "flashing" PMTs that emit light from internal corona 
discharges. Flasher events were removed by two different 
methods. (1) Typical flasher events have broader PMT tim- 
ing distributions than the neutrino events. Events with broader 
timing distributions were eliminated. (2) Since flasher events 
have a tendency to be repeated with similar spatial hit distri- 
bution, the pattern information of observed charge was used 
to eliminate these events. A correlation parameter based on 
the charge pattern was calculated with other data events and 
a "matched" tag was assigned for highly correlated events. A 
cut was applied based on maximum correlation value and the 
number of "matched" with other events. 

(v) Two further event types are eliminated in the fifth re- 
duction step. (1) Events are removed which have > 10 OD 
hits in 200 nsec coincidence preceding the trigger time ( 
—8900 ^ —100 nsec ); this eliminates decay electrons from 
invisible cosmic ray muons that are below Cherenkov thresh- 
old in the inner detector. (2) Cosmic ray muons are removed 
using a more precise fitter and the same criteria as (1) of (iii). 

(vi) Finally, the vertex was required to be within a fiducial 
volume, 2 meters from the wall of the inner detector, and the 
visible energy was required to be greater than 30 MeV. 

TableUshows the number of events for each reduction step. 
Also shown are the the number of Monte Carlo events for each 
reduction step. 



The Super-Kamiokande event sample consists mainly of 
downward-going cosmic ray muons and low energy radioac- 
tivity from parents such as radon. Owing to the double struc- 
ture of the inner and outer detectors, cosmic ray muons are 
easily removed with high efficiency by requiring little or no 
activity in the outer detector For atmospheric neutrino anal- 
ysis, we then consider only events with visible energy above 
30 MeV, where visible energy (Eyis) is defined as the energy 
of an electromagnetic shower that gives a certain amount of 



2. Event Reconstruction 

The fully-contained events underwent a series of recon- 
struction steps in order to classify their origin and properties. 
First, the vertex position of an event was determined using 
PMT hit times; the point which best fit the distribution of PMT 
times (when adjusted for the time of flight of the Cherenkov 
light) was defined as the vertex position. This vertex was re- 
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Reduction step 


Data 


Monte Carlo 


Trigger 


1889599293 


14013.9(100.00%) 


First reduction 


4591659 


14006.3 (99.95 %) 


Second reduction 


301791 


14006.1 (99.94%) 


Third reduction 


66810 


13993.3 (99.85 %) 


Fourth reduction 


26937 


13898.1 (99.17%) 


Fifth reduction 


23984 


13895.3 (99.15%) 


Fiducial volume and 


12180 


13676.7 (97.59%) 


visible energy cuts 







TABLE I: Number of events after each reduction for fully-contained 
events during 1489 days of the detector live- time. The Monte Carlo 
numbers and efficiencies down to the fifth reduction are for events 
whose real vertex is in the fiducial volume, the number of outer de- 
tector hits fewer than 10 and the visible energy larger than 30 MeV. In 
the last line, the fitted vertex is used for both data and Monte Carlo. 



constructed again after particle identification was established, 
to correct for particle track length. The vertex resolution was 
estimated to be 30 cm for single-ring fully-contained events. 
The distribution of vertex position for both data and MC as 
a function of the z-coordinate and -coordinate are shown 
in Figures^] (a) and (b), respectively. In these and several 
further figures, the original Monte Carlo prediction is modi- 
fied by the oscillation of CC interactions according to the 
best-fit parameters (sin^29 = 1.0, Am^ = 2.1 x lO^^eV^), as 
found in Section VI. Only a simple survival probability sup- 
pression is applied for these comparisons, the adjusted sys- 
tematic terms that will be described in Section lVTI 

After an initial ring direction and vertex were found by use 
of the timing method, a Hough transform |76| based tech- 
nique was applied to automatically determine the number of 
Cherenkov rings in an event and their directions. The tech- 
nique was iterative. A second ring was searched for by choos- 
ing possible ring directions based on the Hough map, and a 
likelihood technique was used to determine if a second ring 
from this list of possible rings was more consistent with the 
data than just one ring. If a second ring was found to be nec- 
essary, then this procedure was repeated as often as needed 
(to a maximum of 5 found rings), each time fixing the previ- 
ously found rings, until finally no further rings were necessary 
to fit the data. Figure shows the likelihood difference be- 
tween the 2-ring assumption and a 1-ring assumption. A cut 
was made at likelihood difference of to separate single and 
multi-ring events. The likelihood distributions, especially the 
one for multi-GeV energy region, have a slight difference in 
the peak positions between the data and the Monte Carlo. This 
difference is taken as a source of the systematic error in the 
measurements of the v^, and Vg rates. More details will be dis- 
cussed later Figure^] shows the distribution of the number 
of reconstructed Cherenkov rings for both the data and MC. 

The efficiency for identifying charged current (CC) quasi- 
elastic VeiVfj) events as single-ring was 93.2 (95.8) %, and the 
angular resolution for these single-ring events was estimated 
to be 3.0° and 1.8° for single-ring e-like and /j-like events. 
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FIG. 17: The distribution of fully-contained event vertices in the (a) 
z-coordinate and (b) r^-coordinate, comparing SK data and atmo- 
spheric neutrino Monte Carlo. The points show the data and the his- 
togram shows the Monte Carlo prediction. The Monte Carlo includes 
neutrino oscillation with (sin^28 = 1.00, Am^ = 2.1 x 10"-' eV^). 



respectively. 

To determine the identity of the final state particles, a 
particle identification algorithm was applied which exploited 
systematic differences in the shape and the opening an- 
gle of Cherenkov rings produced by electrons and muons. 
Cherenkov rings from electromagnetic cascades exhibit a 
more diffuse light distribution than those from muons. Fig- 
ures |20| and |^ show observed single-ring e-like and /j-like 
events, respectively. The opening angle of the Cherenkov 
cone, which depends on p(= v/c), was also used to sep- 
arate electrons and muons at low momenta. The validity 
of the method was confirmed by a beam test experiment at 
KEK iIttIi . The misidentification probabilities for single-ring 
muons and electrons were estimated to be 0.7 % and 0.8 % re- 
spectively, using simulated CC quasi-elastic neutrino events. 
The distribution of the likelihood variable used to discrimi- 
nate single-ring electrons and muons are shown for both the 
data and MC for the sub-GeV and muIti-GeV samples in Fig- 
ure|22] In both of these cases there is a clear separation of the 
likelihood variable. 

Figure shows the likelihood variable distribution for the 
brightest ring of FC multi-ring events. Due mostly to overlap- 
ping of Cherenkov photons from multiple particles, the sepa- 
ration of the particle type for a Cherenkov ring in a multi-ring 
event is not as good as that for a single-ring event. 

The identification efficiency was checked using cosmic ray 
muons that stop in the detector and subsequently decay to 
electrons. These events are easily selected by their timing sig- 
nature. The resulting misidentification probabilities for stop- 
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FIG. 18: The distribution of the likelihood difference between a 
single-ring assumption and a multi-ring assumption for sub-GeV 
(top) and multi-GeV (bottom) FC events. The points show the data 
and the histograms show the Monte Carlo prediction. The Monte 
Carlo includes neutrino oscillation with (sin^29 = 1.00, Am^ = 
2.1 X 10~^ eV'^). The hatched histograms show the charged current 
quasi-elastic interactions. 



FIG. 20: An example event display of a single-ring e-like event. Each 
small circle represents a hit PMT and the size of the circle represents 
the number of photons to hit it. In this event the boundary of the 
Cherenkov light is smeared over many PMTs as the light comes from 
numerous positrons and electrons in the electromagnetic shower. 




Number of rings 



FIG. 19: The distribution of the number of identified Cherenkov 
rings, comparing SK data and atmospheric neutrino Monte Carlo. 
The Monte Carlo includes neutrino oscillation with (sin^29 = 1.00, 
Am- = 2. 1 X 10^^ eV^). The hatched histogram shows the charged 
current quasi-elastic interactions. 



ping cosmic ray muons and decay electron light patterns were 
0.4 ± 0.1 % and 1.8 ± 0.5 % respectively, in good agreement 
with the Monte Carlo estimates. This check was performed 
continuously during data-taking, and particle identification 
performance remained stable despite water transparency that 
varied from about 90 m to 120 m. 

Next, the Cherenkov rings were re-fit taking into account 



the expected light pattern given by the particle identifica- 
tion, and in the case of single-ring events, a specialized event 
fitter was applied. After the rings were re-fit and the to- 
tal photo-electrons in the event were apportioned between 
all of the rings, each ring was assigned a momentum based 
on proportion of Cherenkov photons. The momentum of 
a particle was determined from the total number of p.e.'s 
within a 70° half-angle cone relative to the track direction, 
with corrections for light attenuation and PMT angular ac- 
ceptance. The re sulting mom entum resolution is estimated 
to be 0.6 + 2.6/-\/P(GeV/c) % for single-ring electrons and 
1.7-|-0.7/^P(GeV/c)% for single -ring muons. A final pro- 
cedure was performed which utilized the final energy and an- 
gle information of the rings to remove rings which were most 
likely not real. 

Although decay electrons are not used in this oscillation 
analysis, they are a useful signature in other atmospheric neu- 
trino analyses and the search for proton decay; therefore, we 
document their treatment here. Decay electrons were iden- 
tified either as: (a) PMT hits within the same time window 
as the primary event trigger (up to 900 ns later) or (b) a later 
independent event trigger. In the first case, a sliding search 
window of width 30 ns began 100 ns after the primary trigger; 
a decay electron was counted if 40 hits were found in coin- 
cidence above the background level. In the second type, 60 
hits were required in a 50 ns time window, and goodness-of- 
fit for a Cherenkov ring pattern is required. In both cases the 
vertex is known from the primary event and is used to sub- 
tract the time-of-flight of the Cherenkov light. If the decay 
occurs around 900 ns, the hits may be split between the pri- 
mary event trigger and a subsequent event trigger. In some 
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the multi-ring events. Points show the data and the histograms show 
FIG. 21: An example event display of a single-ring ^;-like event. In the Monte Carlo results. The Monte Carlo includes neutrino oscil- 
this event the boundary of the Cherenkov light is sharp as the muon lation with (sin^2G = 1.00, Am^ = 2.1 x 10"^ eV^). The hatched 
travels relatively straight as it comes to a stop. Distant hit PMTs histograms show the charged current interactions, 
come from scattered light and Cherenkov light from delta-rays. 
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FIG. 22: The distribution of particle identification likelihood for 
sub-GeV (top) and multi-GeV (bottom) FC single-ring events, com- 
paring SK data (points) and atmospheric neutrino Monte Carlo 
(histograms). The Monte Carlo includes neutrino oscillation with 
{sin^lQ = 1.00, Am^ = 2.1 x 10"^ eV^). The hatched histograms 
show the charged current interactions. 



analyses, electrons in the time interval 800 ns to 1200 ns after 
the primary trigger are excluded, owing to this splitting effect 
as well as a reduced efficiency due to electrical reflection on 
the PMT cables. The contamination level for these criteria 
is very good, with no events out of 32000 stopping cosmic 
ray muons having more than one decay electron. The effi- 
ciency for fully-contained sub-GeV neutrino interactions was 



3. Background and Efficiency 

The main sources of the background for the FC sample are 
cosmic ray muons, neutrons generated by high energy cos- 
mic ray muons and PMT flasher events. The contamination of 
the background events was estimated for lower energy (visible 
energy lower than 1.33 GeV, called iM^-GeV) and higher en- 
ergy (visible energy higher than 1.33 GeV, called multi-GcY) 
samples separately, since the contamination could have an en- 
ergy dependence. The cosmic ray muon background contam- 
inations to the final FC sample were estimated to be 0.07 % 
for sub-GeV and 0.09 % for multi-GeV jU-like events. They 
were estimated using the distribution of the distance of the 
vertex position from the inner detector wall along the particle 
direction. The neutron background contamination was esti- 
mated to be 0. 1 % for sub-GeV e-like events and multi-GeV 
e-like events. This was estimated using the distribution of the 
distance of the vertex position from the inner detector waU 
(Dwaii)- The contaminations from PMT flasher events were 
estimated to be 0.42 % for sub-GeV e-like events and 0.16 % 
for multi-GeV e-like events using the goodness of the vertex 
fitting and the D„aii distributions. 

The reduction efficiency was estimated using the atmo- 
spheric neutrino MC. The detection efficiency for events 
which satisfy the data reduction conditions (i)-(iv) was 
99.15 % for events which have a true vertex in the fiducial vol- 
ume, Eyis > 30 MeV, and less than 10 hits in the largest OD 
cluster. The systematic error in the event reduction was esti- 
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mated to be 0.2%. The main source of the systematic error 
in the event reduction was the flasher cut based on the pat- 
tern matching algorithm. The systematic error for this cut was 
estimated by mixing the different flasher samples in the atmo- 
spheric MC and comparing the reduction efficiencies. The in- 
efficiency was estimated to be 0.7 % ± 0.2 %, in which 0.2 % 
was considered as the systematic uncertainty. The system- 
atic errors for other reduction steps were negligibly small be- 
cause the reduction efficiency was almost 100% and the dis- 
tributions of cut parameters for data agree with that of MC. 
In the early stage of the experiment, an independent data se- 
lection and reconstruction program was also employed; the 
results were compared and event samples and classifications 
were found to agree by around 95 % LU . 



B. Partially-Contained Events 

1. Data Reduction 

The data reduction for PC events differed from the reduc- 
tion for FC events because of the presence of additional hits in 
the OD. Because these extra hits result from the exiting par- 
ticle (usually a muon), a simple criterion based on the num- 
ber of hit OD tubes could not be used to reject cosmic ray 
background. The criteria used to identify partially-contained 
events are as follows: 

(i) Low energy events with fewer than 1000 total p.e.'s in 
the inner detector were removed, corresponding to muons 
(electrons) with momentum less than 310 (110) MeV/c. By 
definition, an exiting particle in the PC sample must have 
reached the OD from the inner fiducial volume, and so must 
have had a minimum track length of about 2.5 m (correspond- 
ing to muons with > 700 MeV/c momentum). 

(ii) Events for which the width of the time distribution of 
hits in the OD exceeded 260 nsec were rejected, as well as 
events with two or more spatial clusters of OD hits. These cuts 
eliminated many through-going muons, which typically left 
two well separated clusters in the OD. Muons which clipped 
the edges of the detector were eliminated based upon the 
topology of the OD cluster Cosmic ray muons which en- 
tered and stopped in the inner volume of the detector were 
eliminated by excluding events with a relatively small number 
of ID photoelectrons near the OD cluster (1000 p.e.'s within 
2 m). This cut did not remove PC neutrino events because PC 
events produced large numbers of photoelectrons (typically 
3500 p.e.'s) in the region where the particle exited. 

(iii) In the next step, a simple vertex fit and p.e. weighted 
direction estimation were used. A requirement of < 10 hits in 
the OD within 8 m of the back-projected entrance point was 
imposed. Also in this step, flasher events were removed by 
using their broader timing distribution feature. 

(iv) The remaining background still had muons which left 
few or no entrance hits in the OD. These events were rejected 
by requiring the angle subtended by the earliest inner detector 
PMT hit, the vertex, and the back-projected entrance point 
be > 37°. Remaining corner clipping muons were rejected by 
requiring a fitted vertex at least 1.5 m away from the corners of 



the ID volume. A through-going muon fitter was also applied 
to reject events with a well fitted muon track greater than 30 m 
long. 

(v) In the last reduction step, various remaining background 
events were eliminated by several selection criteria: (1) Fully- 
contained events were eliminated by requiring that PC events 
have more than 9 hits in the most highly charged cluster in 
OD; (2) A minimum requirement of 3000 total p.e.'s in the 
inner detector, which corresponds to 350 MeV of visible en- 
ergy, well below that of any exiting muon, was applied to get 
rid of low energy background events; (3) Clusters in the OD 
were searched for again with the same clustering algorithm 
used in the 2nd reduction step but with different clustering pa- 
rameters. Events were eliminated if there existed two or more 
clusters with more than 10 p.e.'s and they were apart by more 
than 20 meters. Some obvious through-going muons were re- 
moved by this cut; (4) After those steps, most remaining back- 
ground events are due to the imprecision of the fast fitters used 
to quickly filter the data stream. A precise fitting algorithm 
was then applied to obtain more accurate information on ring 
direction and vertex position. With much more accurate infor- 
mation of the event, we were able to eliminate most remaining 
through-going and stopping events based on their distinct ge- 
ometry and OD signatures. (5) Some through-going muons 
have a very special geometry — they passed through the tank 
vertically along the wall of ID. These events were eliminated 
by counting the number of p.e.'s and hits in the OD within 
the region defined by an 8 m sphere around the top and bot- 
tom fringes and checking the time interval between the av- 
erage timings of those top and bottom hits. (6) Remaining 
cosmic ray muon background events are those entering the 
ID through relatively weak OD regions — there are four holes 
covered by veto counters on the top of OD through which ca- 
bles run. Events with a veto counter hit were eliminated, as 
well as those satisfying a detailed cosmic ray muon consis- 
tency requirement. 

(vi) After this final reduction step, events were scanned by 
physicists to check the data quality. However, no event was re- 
jected based on the scanning. Finally, the vertex was required 
to be within a fiducial volume, 2 meters from the wall of the 
inner detector. The final event sample is an almost 100 % pure 
V sample. The background contamination has been estimated 
to be about 0.2 %. 

Table HH shows the number of events after each reduction 
step and the detection efficiency of PC events as a function of 
reduction steps. 



2. Event Reconstruction 

The partially-contained events were reconstructed using in- 
ner detector PMT information by similar vertex, direction fit 
and ring-counting algorithms, as were applied to fully con- 
tained events. 

For some PC events, however, the direction fit was slightly 
modified. Under some conditions, rather than using the re- 
sults of the precise fitting algorithm for the PC event direction, 
outer detector spatial information was used instead. In order 
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Reduction step 


Data 


Monte Carlo 


Trigger 


1889599293 


1,417.0(100.0%) 


First reduction 


34536269 


1,402.8 (99.0%) 


Second reduction 


5257443 


1,334.7 (94.2%) 


Third reduction 


380053 


1,318.7 (93.1 %) 


Fourth reduction 


53825 


1,246.2 (87.9%) 


Fifth reduction 


1483 


1,201.0 (84.8%) 


Fiducial volume 


911 


1,129.6 (79.7%) 



TABLE 11: Number of events after each reduction step for partially- 
contained events during 1489 days of the detector live-time. The 
Monte Carlo efficiencies are for events whose real vertices are in the 
fiducial volume and the number of outer detector hits more than 9. 
In the last line, we used the events whose fitted vertices are inside 
the fiducial volume both for data and Monte Carlo. The Monte Carlo 
does not include neutrino oscillation. 
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FIG. 24: The distribution of partially-contained event vertices in the 
(a) z-coordinate and (b) r^-coordinate, comparing SK data and atmo- 
spheric neutrino Monte Carlo. The Monte Carlo includes neutrino 
oscillation with (sin2 2e = 1.00, Am^ = 2.1 x 10"^ eV^). 



to use the outer detector cluster for the direction, the number 
of tubes in the largest OD cluster was required to be greater or 
equal to 20. If this condition was satisfied, and also, if the ID 
PMT nearest to the projected ID exit point of the fitted track 
with more than 200 p.e.'s was more than 200 cm away (i.e. no 
clear ID exit point), or if there was a clear exit point in the ID 
but yet there were more than 800 ID PMT hits with more than 
200 p.e.'s in each of them (i.e. saturated in our electronics), 
then the vector from the fitted vertex to the largest OD cluster 
was used for the PC event direction. Otherwise, the standard 
direction provided by the precise fitting algorithm was used. 



The estimated vertex position resolution for PC events was 
64 cm. The angular resolution for the penetrating particle in a 
PC event was estimated to be 2.8°. 

Finally, the fiducial volume cut was applied. The event rate 
in the fiducial volume was 0.62 events/day. 



3. Background and Efficiency 

The background for the PC sample originates from cosmic 
ray muons. They were efficiently removed by the reduction 
steps mentioned above. We estimated the contamination of 
non-neutrino events in the fiducial volume by two methods. 
One method utilized scanned results. After applying all the 
reduction steps, all the events were scanned and the estimated 
fraction of background contamination was found to be 0.2 % 
in the fiducial volume. Another method was by examining 
the vertex distribution of non-neutrino events as a function of 
distance from the wall. By extrapolating the distribution from 
outside the fiducial volume, we obtained that the contamina- 
tion of the background was less than 0.1 %. Since these two 
results were statistically consistent, we took the larger num- 
ber (0.2 %) as the contamination of background events in the 
fiducial volume. Figures |24] (a) and (b) show the distribution 
of the vertex position for both data and MC as a function of 
the z-coordinate and r^-coordinate. Some contamination of 
background is evident near the side and top PMT walls. How- 
ever, no evidence for substantial background contamination is 
seen in the fiducial volume. 

The PC reduction efficiency was estimated based on Monte 
Carlo events. The definition of partially-contained events is 
that ( 1 ) the interaction point of the parent neutrino is inside the 
fiducial volume, and (2) number of outer detector hits within 
8 m around an estimated exiting point is larger than nine. We 
applied those five reduction steps to the atmospheric neutrino 
Monte Carlo sample and after each reduction step, we counted 
the number of PC events left inside fiducial volume and cal- 
culated the efficiency of this particular reduction step. We 
found the overall the efficiency of the reduction for PC neu- 
trino events was 79.7 %. 

The systematic uncertainties on the reduction efficiency 
from the first to the fourth steps were estimated by two meth- 
ods. For the ID contribution, they were estimated by compar- 
ing the distributions of cut parameters used in the reduction 
criteria for the data and Monte Carlo. For the OD part, we 
created Monte Carlo samples with two different sets of OD- 
related tuning parameters. Then, the change in reduction ef- 
ficiency is the OD-related systematic error. The two sets pro- 
duced different amounts of light in the OD within the limits 
of good overall agreement with the standard tuning sample 
(several hundreds of well-measured stopping muons). The 
estimated uncertainties were 1.5% and 1.4% for inner- and 
outer-detector related selection criteria, respectively. 

The systematic uncertainty in the fifth reduction efficiency 
is mainly from ID variables and involves the precise fitting 
algorithm. The main contributions come from the cluster cut 
mentioned above. The uncertainty on the cluster cut was es- 
timated by examining the distributions of cut variables. The 
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uncertainties in other cuts were assumed to be the inefficiency 
in each cut, since the inefficiency itself was small compared 
with other errors. Combining these estimated uncertainties, 
we obtained the value of 1 .6 % for the fifth reduction step. 

Thus the overall systematic uncertainty in the reduction of 
the PC events was estimated to be 2.6 %. 



C. Upward-going muon Data 

1. Data Reduction 

The upward-going muons observed in Super-K are classi- 
fied into two categories: (1) upward stopping muon events 
having only an entrance signal in the OD; and (2) upward 
through-going muon events having both entrance and exit sig- 
nals in the OD. The criteria used in Sec. lVB Ti to determine an 
event entry or exit were used: 10 OD hits in-time and within 
8 m of the tracks projected entry or exit point constitutes a 
muon entry or exit signal. We required the geometrical tra- 
jectory of through-going muons to be greater than 7 m in the 
inner detector, and we imposed an equivalent 7 m path-length 
cut on upward stopping muons based on the momentum re- 
construction using Cherenkov light. 

The purpose of the data reduction is to isolate the upward 
muon events and the horizontal muon events (needed for back- 
ground estimation), provide a classification of stopping or 
through-going muon type, and to reject the background from 
cosmic ray muons and noise such as flashing PMTs. Decay 
electrons associated with stopping muons were also saved. A 
charge cut of 8,000 < g < 1,750, 000 p.e.'s in the ID was ap- 
plied. For a muon, Q{p.e.) ~ 25 ■ L{cm), ensuring that that we 
can detect all muons with path-length > 7 m while eliminating 
events at lower energies. At very high ID charge correspond- 
ing to ~ 1 , 750, 000 p.e.'s the ID electronics becomes saturated 
causing the muon fitters fail. 

To isolate the rate of about one neutrino induced upward- 
going muon per day from the remaining background of about 
2x10^ cosmic ray muons, we used a logic tree involving as 
many as seven different muon fitters. Some of these fitters 
were specialized to fit stopping muons, others were special- 
ized for through-going muon events, and some of them were 
specialized to fit muon events with Bremsstrahlung. The main 
logic behind the upward-going muon reduction was that if a 
muon fitter classifies an event as upward with a goodness of 
fit which was above the fitter's goodness threshold then the 
event was automatically saved. Conversely, if a muon fitter 
classified an event as downward with a goodness of fit which 
was above the threshold then the event was automatically re- 
jected. If a muon fitter classified an event as horizontal and 
with goodness above threshold, or if the fitter could not give 
a good fit for the event, the event was passed to the next fit- 
ter. This sequence continued until all the events had passed 
through all the fitters or had been classified. If no fitter was 
able to give a good fit then this event was automatically re- 
jected. If at least one fitter classified this event as horizontal 
then the event was saved. All events from the output of the 
upward muon reduction were then passed to the precise fitter 



which is described in Sec. lVC2l 



2. Event Reconstruction 

All events from the output of the upward-going muon re- 
duction were passed through the precise fitter. The basic algo- 
rithm was identical to that used for the vertex and direction 
determination for single ring fully-contained and partially- 
contained events. The fitter assumes that the particle is a 
muon and the vertex position of the event is at the inner de- 
tector surface. However, when the muon produces an ener- 
getic electro-magnetic shower, the assumption of single non- 
showering muon does not give an accurate direction. For these 
events, the information of OD hit is used to determine the par- 
ticle direction. The angular resolution of the fitter was about 
1.0° for both through-going and stopping muons. Taking into 
account multiple scattering from the point of muon creation, 
68 % of through-going and stopping muons fit within 1.3° and 
2.4° of the muon's true initial direction, respectively. The di- 
rection determined by this fitter was used in the neutrino os- 
cillation analysis. 



3. Background and Efficiency 

The effective detection efficiency for the data reduction pro- 
cess was estimated by a Monte Carlo simulation, and was 
found to be 102.4 % for upward stopping muons and 95.9 % 
for upward through-going muons. The efficiency higher than 
100% for upward stopping muons is due to a slight bias in 
the separation of stopping and through-going muons, causing 
a small fraction of the more numerous through-going muons 
to be misidentified as stopping muons. The up/down sym- 
metry of the detector geometry allows a check of this Monte 
Carlo efficiency calculation using real cosmic ray induced 
downward-going muons. The efficiency is approximately 
constant for — 1 < cos © < 0; bin-by-bin efficiencies are listed 
in Table lXlfl in the Appendix. 

After the reduction, the precise fitter described in Sec. IVC2l 
was applied to determine the entry position and the muon di- 
rection. This is the final direction used for further physics 
analysis, including the determination of upward versus down- 
ward classification. Events selected as upward by this direc- 
tion were then scanned by eye with a graphical event display 
program in order to reject difficult-to-remove instances of cor- 
ner clipping or bremsstrahlung cosmic ray muons and noise 
events. The event scanning only rejected events judged to be 
background and did not change the direction and the vertex 
decided by the precise fitter, nor the stop/through judgment 
made by the reduction programs using the entrance and exit 
points of this fit. The event scanning was done independently 
by two physicists and testing had shown that both scanners 
had never rejected the same good upward-going muon event. 
About 50 % of the events remaining after all automated reduc- 
tion steps were rejected by this final scan. Table [111] summa- 
rizes the data reduction for upward muons. 
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Reduction step 


Data 


Monte Carlo 




stopping through-going 


stopping through-going 


Trigger 


2129729843 


697.1 (100%) 1741.0(100%) 


Reduction 


89911 


693.9 (99.5 %) 1722.3 (98.9 %) 


Precise fitter (cos 9 < 0) 


4266 


692.4 (99.3 %) 1721.7 (98.9 %) 


Scan 


2447 




Stop-through separation and > 1.6GeV 


458 1856 


713.5(102.4%) 1669.5 (95.9%) 


CR BG subtraction (subtracted ev) 


417.7 (40.3) 1841.6 (14.4) 





TABLE III: Number of events after each step of the data selection for upward muons during 1646 days of the detector live-time. The Monte 
Carlo efficiencies are relative to the generated events with track length longer than 7 m (for through-going muons) or with energy higher than 
1.6 GeV/c (for stopping muons). The efficiencies are for Monte Carlo upward muon events with track length in the inner detector longer than 
7 m (for through-going muons) or with muon momentum at the wall of the inner detector higher than 1.6 GeV/c (for stopping muons). The 
Monte Carlo does not include neutrino oscillation. 



Near the horizon, horizontal cosmic ray muons are a non- 
negligible source of background for both through-going and 
stopping upward muons. Because of finite fitter resolution 
and multiple Coulomb scattering of muons in the nearby rock, 
some downward going cosmic ray muons may appear to be 
coming from cos0 < 0. Figure l25l shows the zenith versus 
azimuth directions for the upward-going muon sample. Clus- 
ters of cosmic ray downward muons are seen for relatively 
thin overburden directions. Figure shows the zenith an- 
gle distribution of upward muon candidates near the horizon 
for two different regions in azimuth. The thick overburden 
region has negligible downward going cosmic ray muon con- 
tamination, even above the horizon. The thin overburden re- 
gion has non-negligible contamination. The shape of the dis- 
tribution above the horizon was extrapolated below the hori- 
zon to estimate the background contamination in the upward 
muon sample. The number of background events, based on 
Fig-Hi] to the upward stopping muon signal were estimated to 
be UA±7.2{stat)±6.0{sys)md40.3±13.7{stat)±4:3{sys) 
events for the through-going and stopping muon samples, re- 
spectively. Horizontal muon background was contained in the 
—0.1 < cos© < zenith angle bin, and was subtracted from 
this bin. The stopping muon contamination is larger than the 
through-going contamination since lower muon energies al- 
low larger scattering angles. 

Away from the horizon, a potential source of upward parti- 
cles is the photo-production of pions by energetic muons that 
pass nearby the detector |78|. The pathlength requirement of 
7 m within the inner detector limits the background from this 
source to ~ 0.01 % in the upward through-going muon data 
sample and ~ 0.3 % in the upward stopping muon data sam- 
ple. 



4. Expected Upward-Going Muon Signal 

The expected upward-going muon flux was calculated us- 
ing the same tools as for the contained vertex events, extended 
to higher energies and to outside the detector volume. The 
input neutrino flux (see Sec. HI} in Ref. 12 8.1 was used up to 
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FIG. 25: The zenith versus azimuth directions for a 641.4 day sample 
of upward muons used to estimate background contamination. The 
dense regions of events near the horizon correspond to thin regions of 
the mountainous overburden. The upper panel is for upward stopping 
muons and the lower panel is for upward through-going muons. 



1 TeV. At 1 Te V, the calculated flux in Ref. 17911 was rescaled 
to that in Ref. i28il . Above 1 TeV, the rescaled flux in f79il was 
used up to 100 TeV. The target volume for these neutrinos is 
primarily the rock around the detector, parameterized as stan- 
dard rock, with Z=l 1, A=22 and density = 2.7 gm/cm^. How- 
ever, neutrinos interacting in the water of the OD and insensi- 
tive region can also be seen as upward-going muons (1.8 % of 
through-going and 6.6 % of stopping muons), so water inter- 
actions were also simulated. 

The neutrino interactions were modeled as discussed in 
Sec. lIVI The same GEANT detector simulation discussed pre- 
viously was used to track muons from the interaction vertex 
through the rock into the detector itself. The output of the de- 
tector simulation was passed through the same reduction and 
fitting routines as was the real upward-going muon data. 100 
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FIG. 26: The zenith angle distribution of upward muon candidates 
near the horizon for two different regions in azimuth. Regions (1) 
and (2) are thin and thick overburden respectively. The live-time 
of upward and downward events is 1645.9 days and 641.4 days re- 
spectively. The number of downward events is normalized to be 
1645.9 days live-time. The upper panel is for upward stopping 
muons and the lower panel is for upward through-going muons. 



years equivalent exposure was generated. The results of this 
Monte Carlo data before and after the reduction were used to 
estimate efficiencies and systematic errors as well as provide 
an expected upward-going muon signal. TablellUlsummarizes 
the detection efficiency at each step of the data reduction. 



D. Observations 

1. Contained Vertex Events 

We have accumulated 1489.2 days of FC and PC data from 
May 17, 1996 to July 16, 2001. These are compared with sta- 
tistically larger samples of Monte Carlo events based on the 
two neutrino interaction models |47, 48], both equivalent to 
100 years of the detector exposure. However, upward-going 
Monte Carlo muons were only generated based on NEUT. 
Both Monte Carlo samples were generated based on the flux 
model of Ref. |28|. The Monte Carlo samples were processed 
by the same event-selection and event reconstruction steps 
as the real events. Throughout this paper, unless otherwise 
indicated, we used Monte Carlo events generated based on 
the flux model of Ref L28iJ and the NEUT neutrino interaction 
model 1 47]. 

Fully-contained events were divided into two sub-samples 
according to the reconstructed visible energy. We refer to 
the event sample below 1.33 GeV as sub-GeV, and above 
1.33 GeV as multi-GeV. Fully-contained events were further 
divided into the events with single reconstructed Cherenkov 
ring, single-ring, and events with more than one rings, multi- 
ring. All single-ring events were classified as either e-like 



or /7-like based on the PID result. Lower energy cuts were 
applied only to the single-ring sample, Pe > lOOMeV/c for 
e-like and > 200MeV/c for jU-like. In addition, multi- 
ring events were used to study the atmospheric neutrino flux. 
A simple set of criteria that the most energetic ring in a 
multi-ring event was /j-like with P^, >600 MeV/c and Eyjs > 
600 MeV selected relatively pure CC events. The estimated 
fraction of CC in this sample was 90.5 % and 94.9 % for the 
sub- and multi-GeV energy ranges, respectively. On the other 
hand, a similar set of criteria for e-like events obtained only 
a 54.4 % pure CC Vg sample, and therefore we decided to use 
multi-ring /j-like events only. 

Table llVl summ arizes the number of observed events in the 
sub-GeV and multi-GeV samples as well as the expected num- 
ber of events in the absence of neutrino oscillations. The frac- 
tion of various neutrino interaction modes, predicted by the 
Monte Carlo sample, are also listed. Figure |^ shows the 
event rates for contained events as a function of the elapsed 
days. The event rate should change due to the solar modula- 
tion. However, the expected decrease in the event rate from 
minimum solar activity (which approximately corresponds to 
the period when the SK-1 started taking data) to maximum so- 
lar activity (which approximately corresponds to the period 
when the SK-1 finished taking data) period is 6-7 % for sub- 
GeV events and 3-4 % for multi-GeV events. The data cannot 
distinguish a constant event rate from the expected rate change 
due to solar activity. 
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FIG. 27: Event rates as a function of elapsed SK running time for 
contained events. 



Figures [T7| |24| show the reconstructed vertex distributions 



for FC and PC events projected on the z and r 



axes. The shape of the vertex distributions of the data and MC 
agree well in the fiducial volume. Near the fiducial volume 
boundary, the event rate for FC events slightly decreases, and 
the event rate for PC events increases. This is because ener- 
getic muons produced by neutrino interactions near the wall 
tend to escape and are identified as PC events. Peaks near 
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(NEUT) (NUANCE) 






(Flux A) 






(FluxB) (FluxC) 


(Flux A) 


sub-GeV 


8941 


9884.3 






9967.8 


10619.4 


9074.2 


single-ring 


6580 


7092.6 






7273.2 


7643.3 


6694.0 


e-like 


3353 


2879.8 


2533.9 (88.0%) 


66.3 (2.3%) 279.6 (9.7%) 


2944.2 


3069.5 


2762.3 


/j-like 


3227 


4212.8 


22.8 (0.5%) 3979.7 (94.5%) 210.3 (5.0%) 


4329.0 


4573.9 


3931.6 


multi-ring 


2361 


2791.7 






2694.6 


2976.0 


2380.2 


jU-like 


208 


322.6 


11.6 (3.6%) 


292.0 (90.5%) 18.9(5.9%) 


301.5 


342.1 


274.0 


R = 


0.658 ± 0.016 (stat) ± 0.035 (sys) 




0.655 


0.646 


0.676 


multi-GeV 


2901 


3472.0 






3212.6 


3708.7 


3179.3 


single-ring 


1397 


1580.4 






1456.8 


1676.6 


1463.7 


e-like 


746 


680.5 


562.2 (82.6%) 


47.6 (7.0%) 70.7 (10.4%) 


635.3 


729.2 


635.3 


iU-like 


651 


899.9 


3.6 (0.4%) 


894.2(99.4%) 2.1(0.2%) 


821.4 


947.4 


828.4 


multi-ring 


1504 


1891.6 






1755.9 


2032.1 


1715.5 


/j-like 


439 


711.9 


16.6 (2.3%) 


675.8 (95.0%) 19.4 (2.7%) 


645.9 


749.1 


618.9 


partially-contained 


911 


1129.6 


20.8(1.8%) 


1098.8 (97.3%) 10.0 (0.9%) 


1065.0 


1236.6 


1074.9 


Rpc+PC = 


0.702 


+0.032 / . ^ 
-0.030 W'™-' 


± 0.101 (sys) 




0.705 


0.699 


0.699 



TABLE IV: Summary of the sub-GeV, multi-GeV and PC event samples compared with the Monte Carlo prediction based on the neutrino 
interaction model of Ref. l47ll (NEUT) and neutrino flux calculation of Ref. 1 28], as well as different flux models. Fluxes A, B and C refers to 
l23l . and |29], respectively. The Monte Carlo prediction with NUANCE |48] and Flux A is also shown. The Monte Carlo predictions do 
not include neutrino oscillations. 



the edge of the ID for both the MC and the data are caused 
by a constraint of the vertex reconstruction programs: the re- 
construction of the vertex is restricted within the ID, and the 
events whose vertex is estimated outside of the ID are con- 
strained to be within the ID, where they pile up at the edge. 
The peak at z = 1810.0 cm in the distribution for FC data is 
caused by the cosmic ray muons passing through inefficient 
regions of the OD. These muons are safely rejected by the 
fiducial volume cut. 

Fig. shows the reconstructed momentum distributions 
for FC single-ring events. The data and MC show good 
agreement except for significantly fewer numbers of FC /j- 
like events. Fig. 1291 shows the visible energy distribution for 
FC multi-ring /j-like and PC events. The PC data have more 
events than the Monte Carlo prediction at the highest ener- 
gies. This could indicate that the neutrino energy spectrum in 
the Monte Carlo is too soft around 100 GeV. We note that the 
upward through-going data also suggest that the neutrino en- 
ergy spectrum is too soft (see Table lVIl ) The neutrino energy 
spectrum up to 1 TeV was considered in the Monte Carlo pre- 
diction for the FC and PC samples. We estimated that two PC 
events are expected with visible energy above 100 GeV from 
neutrinos above 1 TeV. 

The flavor ratio of the atmospheric neutrino flux, (v^ + 
V^)/(V£. + Ve), is predicted with 3 % accuracy. As shown in 
Table II Vl the particle identification for FC single-ring events 
gives a good estimation of the flavor of the parent neutri- 
nos, and the ratio of the number of e-like events and /j-like 
events, (/j/e), gives a good estimation of the flavor ratio of 
the atmospheric neutrinos. We define a (fj/e) double ra- 



tio, R = {^i/ e)oATA/ {^^1 ^)mc- Without neutrino oscillation, R 
should be consistent with unity. R is measured to be: 

^™fc-G«y =0.658 ±0.016 ±0.035, (2) 

for the sub-GeV sample. 

A substantial fraction of muons in the multi-GeV energy 
range exit from the inner detector and are detected as PC 
events. The partially-contained event sample is estimated to 
be 97 % pure CC V;j interactions, even without requiring any 
particle identification or ring-number cuts. Therefore, we add 
FC single-ring and PC event totals when calculating R in the 
multi-GeV range. We measured R in the muIti-GeV energy 
range to be: 

Rmulti-GeV+PC = ^■^^'^^^1^ ±0.101. (3) 

Systematic uncertainties in the double ratio R have been dis- 
cussed in detail in Refs. fl]] and 101 ■ These eiTors have been 
re-evaluated and are estimated to be 5.3 % for sub-GeV and 
14.4% for muIti-GeV events. The sources of the systematic 
uncertainties in R are listed in Table [V] which include both 
theoretical and experimental errors. Among the experimen- 
tal systematic errors, the separation of single- and multi-ring 
events is the largest source of the systematic uncertainties. As 
shown in Fig. ^] the distributions of the likelihood differ- 
ence between the single-ring and multi-ring assumptions have 
slight shifts in the peak positions between the data and the 
Monte Carlo. These differences could cause systematic un- 
certainties in the number of identified single-ring e-like and 
/J-like events (which are summarized in Table IXt. Since the 



20 




1 r 

! u I I I I I I I I I I I I 

10"^ 1 10 10^ 

Momentum (GeV/c) 

FIG. 28: Momentum distribution of FC single ring e-lilce (top) and /j- 
like (bottom) events. The sharp cut in the muon momentum spectrum 
in the high energy end is due to the requirement on fully containment. 
The points show the data, solid lines show the Monte Carlo predic- 
tion without neutrino oscillation and dashed lines show the oscillated 
Monte Carlo events with {sin^lQ = 1.00, Am^ = 2.1 x 10"^ eV^). 



magnitude of the uncertainty is different between e-like and /j- 
like events, and since the separation of single- and muhi-ring 
is not apphed for PC events, the uncertainty in the single- and 
multi-ring separation causes the uncertainties in the R mea- 
surements. 

Figure l30l shows the expected {fj/e)Data/ {fj/e)MC in the 
presence of neutrino oscillation for sub- and multi-GeV sam- 
ples as a function of Arn^. Data are consistently explained by 
neutrino oscillations for Arn^ in the range of 10-^ to lO-^eV^. 

Figurel^shows the measured and expected numbers of FC 
and PC events as a function of the cosine of the zenith an- 
gle (cos 0), cos = — 1 refers to upward-going and cos 0=1 
refers to downward-going. Single-ring sub-GeV events are 
separately shown in two parts, pi < 400 MeV/c and pi > 400 
MeV/c, where pi is the lepton momentum. In the momen- 
tum range below 400 MeV/c, the angular correlation between 
the neutrino and outgoing lepton is very poor, the shape of 
the atmospheric neutrino flux is largely washed out, and the 
zenith angle distributions for the charged leptons should be 
approximately flat. Fig. |32] shows the angular resolution of 
the neutrino directions as a function of the momentum. The 
angular resolution is defined as the angular difference between 
the parent neutrinos and the reconstructed directions in which 
68% of the events are included. 

We have also studied the azimuthal dependence of the at- 
mospheric neutrino data. This is a sensible consistency check. 
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FIG. 29: Visible energy distribution of multi-ring /j-like events (top) 
and PC events (bottom). The points show the data, solid lines show 
the Monte Carlo prediction without neutrino oscillation and dashed 
lines show the oscillated Monte Carlo events with (sin^29 = 1.00, 
Am^ = 2A X 10"^ eV^). 



as neutrino oscillation should not cause any azimuthally de- 
pendent deficit since all path lengths at a given zenith an- 
gle are equal. In a well-selected data sample, the azimuth 
rates exhibit the famous east-west effect, which was used in 
the 1930's to demonstrate that cosmic rays were positively 
charged. The effect is caused by the deflection of primary pro- 
tons in the earth's magnetic field, where trajectories from the 
east are blocked by the bulk of the earth. This also results in 
a deficit of atmospheric neutrinos arriving from the east, with 
the strongest effect at the lowest neutrino energy. We max- 
imized our sensitivity to this effect by selecting lepton mo- 
menta between 400 and 3000 Me V/c in the zenith angle range 
cos0| < 0.5. The low momentum cut ensures good pointing 
resolution, and the high momentum cut diminishes the con- 
tribution from high energy primary protons that are insuffi- 
ciently deflected. The zenith requirement enhances the statis- 
tical sensitivity as the depletion only occurs near the horizon. 
Figure |33l updates our previous result llsoll to the final data re- 
ported here. That our data reproduces the prediction for this 
effect implies that the model for geomagnetic cutoff in the 
flux prediction is accurately accounted for, and checks the ba- 
sic features of neutrino production and scattering. 

While the measured shape of zenith angle distributions for 
e-like events is consistent with expectations, both the FC 
/j-like and PC samples exhibit significant zenith-angle de- 
pendent deficits. The up-down ratio, where U is the num- 
ber of upward-going events (— 1 < cos0 < —0.2) and D is 
the number of downward-going events (0.2 < cos0 < 1), 
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"Different flux calculations predict different energy dependences that cannot 
be explained by a simple spectral index uncertainty. See lower part of Fig. Ill 
Uncertainty of the relative normalization of the fully-contained multi-GeV 
and partially-contained sample gives the systematic error in double ratio R. 

''The mean free path of hadrons in was changed by 30%. Also the 
uncertainty in the pion energy spectrum produced by neutrino interactions, 
defined to be the difference between the interaction models A and B, is taken 
into account. 

TABLE V: Sources of the systematic errors in double ratio R{= 
(pi e)Dfj:ji^l (}xl e)Mc) for the sub-GeV and multi-GeV samples. Esti- 
mated uncertainty in each source of the systematic error is described 
in Tables lVnIlvnilliXl andlXl 



is measured to be: U jD = 1.133+JJ JJg^^ ± 0.009 for single- 
ring sub-GeV e-like events in the momentum range below 
400 MeV/c, V ID= 1.082+°°^ ± 0.024 above 400 MeV/c, 

' — U.UoU 

U ID = 0.964+°°^^ _^ 0.008 for single-ring sub-GeV ^-like 
events below 400 MeV/c, = 0.670+ ° °3^^ ±0.012 above 
400 MeV/c, UlD = 0.961+°°7^^ ±0.016 for the multi-GeV e- 
like events and U lD = 0.55 1+° °3^3^ ±0.004 for the single-ring 
multi-GeV /j-like plus PC events. 

Many systematic uncertainties are canceled for the up- 
down ratio and the remaining sources of the uncertainties are: 
uncertainty in the flux calculation, 0.5 % and 0.8 % for sub- 
GeV e-like and /7-like events in the momentum range below 
400 MeV/c, 2.1% and 1.8% for e-like and ^u-like events above 
400 MeV/c, and 1.5 % and 0.6 % for multi-GeV e-like events 
and multi-GeV /j-like events plus PC events; uncertainty in the 
angular dependence of absolute energy calibration, 0.5 % and 
0.2 % for sub-GeV e-like and /j-like events in the momentum 
range below 400 MeV/c, 0.4 % and 0.4 % for e-like and ^-like 
events above 400 MeV/c, and 0.8 % and 0.4 % for multi-GeV 



FIG. 30: Expected (^1 e)Datal (pl ^^MC for singe-ring sub- and multi- 
GeV -I- PC samples as a function of Am^ for full ^ mixing. The 
values for the data together with ±la statistical errors are shown by 
the horizontal lines. The systematic errors are shown by the band in 
the expectation. 



e-like events and multi-GeV /j-like events plus PC events; and 
uncertainty in the non-neutrino background such as cosmic 
ray muons, <0.4%, <0.1 %, <0.2%, <0.2% for sub-GeV 
e-like, /j-like, multi-GeV e-like, /j-like plus PC events, respec- 
tively. In total, the systematic uncertainty for is 0.8 % 
and 0.8 % for sub-GeV e-like and /j-like events in the mo- 
mentum range below 400 MeV/c, 2.2 % and 1.8 % above 400 
MeV/c, 1.7 % and 0.7 % for multi-GeV e-like and ^-like plus 
PC events, respectively. While the ratio for e-like events is 
consistent with 1, the /j-like up-down ratio for the multi-GeV 
data differs from 1 by more than 12 standard deviations. Fig- 
ure |33 shows the expected V ID ratios as a function of Am^. 
Data are consistently explained by neutrino oscillations with 
Aw^ in the range of 10^^^ to lO^-^eV^. 



2. Upward Muon Events 

The upward-going muon data used in this analysis were 
taken from May 1996 to July 2001. The detector live-time 
was 1645.9 days. Though spanning the same period of calen- 
dar time, this live-time was larger than that of the contained 
vertex events because the reconstruction of long path length 
muons is less sensitive to detector conditions, allowing looser 
run selection criteria. Fig.|35lshows the event rates as a func- 
tion of the elapsed days for upward going muons. The event 
rates for these samples are stable. Table IVll summarizes the 
number of observed events in the upward-going muon data 
sample and the corresponding flux and expected flux. 

The systematic errors on the observed number of events 
compared to the Monte Carlo predictions are: the 2 % en- 
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FIG. 31: The zenith angle distribution for fully-contained 1-ring events, multi-ring events, partially-contained events and upward muons. The 
points show the data, box histograms show the non-oscillated Monte Carlo events and the lines show the best-fit expectations for ^ Vt 
oscillations with sin^29 = 1.00 and Am^ = 2.1 x 10^^ eV^. The best-fit expectation is corrected by the 39 systematic error terms, while the 
correction is not made for the non-oscillated Monte Carlo events. The height of the boxes shows the statistical error of the Monte Carlo. 



ergy scale uncertainty leads to a ^ % error in the stop- 
ping muons due to the 1.6 GeV/c cut; the reduction effi- 
ciency for stopping (through-going) muons has an uncertainty 
°f -1 25 (-0 54 '^^^ stopping/through-going separation 
-0 38 (where "H-" means through-going muons misidenti- 
fied as stopping). As in the contained event analysis, com- 
parison of data and expectations is done between observed 
number of events and the live-time-scaled MC number of 
events. However, to facilitate comparisons with other exper- 



iments, these numbers are also presented in units of flux as 
described in f3!,' 4]. The additional systematic uncertainty in 
the observed through-going (stopping) flux comes from effec- 
tive area of 0.3 % and the live-time calculation (0.1 %). The 
absolute expected flux has theoretical uncertainties of at least 
20 % in the normalization for high energy (> 100 GeV) neu- 
trinos and 5 to 10 % from interaction model differences. 

The zenith angle distributions of the upward through-going 
and stopping muons are shown in Fig. |^ The shape of 
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event class 


# events 


# expected 
(Flux A) 


flux 

(xl0-i\™-2^-l^r-i) 


expected flux 
(xlO-'3c-OT-2^--liT-i) 


# expected 
(Flux B) (Flux C) 


^thru 
%. = ^aopl^thru 
*«o/, +*,/,„, 


417.7 
1841.6 
0.227 
2259.3 


713.5 
1669.5 
0.427 
2382.9 


0.381 ±0.024lJ5;jjJ]^ 
1.661 ±0.040+J5o[^ 
0.229 ±0.015 ±0.003 
2.042 ±0.046l[J;|J[^ 


0.648 ±0.145 
1.506 ±0.337 
0.430 ±0.058 
2. 154 ±0.482 


681.5 790.5 
1644.3 1974.9 

0.414 0.400 
2325.8 2765.4 



TABLE VI: Summary of observed and expected results for upward-going muons during 1645.9 live-days. The first and the second errors in the 
observed flux show statistical and systematic errors, respectively. Expected event rates based on different flux models are also shown. Fluxes 
A, B and C refer [28], [25 J and [22J, respectively. 
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FIG. 32: Angular resolution of the neutrino direction as a function 
of the outgoing charged lepton momentum. The angular resolution is 
defined as the angular difference between the parent neutrino and the 
reconstructed direction for which 68 % of the events are included. 
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FIG. 33: The azimuth distribution of the sample of events selected 
for analysis of the east-west effect. The Monte Carlo histogram is 
normalized to the total number of the real data. 0, 90, 180 and 270 
degree azimuthal angles show particles going to north, west, south 
and east, respectively. 



the zenith angle distribution is predicted accurately. There- 
fore, the vertical to horizontal ratio was taken to study the 
effects of neutrino oscillations, where V and H represent the 
number of through-going events with — 1 < cos® < 0.5 and 
—0.5 < COS& < 0, respectively. The V/H ratio for the data 
was 0.497 ±0.022(sfflf) ±0.003(sys), while the no-oscillation 
prediction was 0.586 ± 0.019(iyi). Taking into account sta- 
tistical, systematic, and theoretical uncertainties (using the 
methods discussed in Section fvH . the V/H ratio of the up- 
ward through-going muon sample was smaller than the pre- 
diction by 3 standard deviations. The observed flux falls off 
much more rapidly than predicted as the zenith angle ap- 
proaches the nadir. Fig.|36lshows the expected and observed 
V/H ratio of the upward through-going muon events. The 
observed ratio suggests that Am^ lies in the range of either 
(1 -3) X 10-^ or (5 - 10) X 10-2 eV^. 

The large uncertainty in the absolute flux normalization can 
be canceled by taking the stopping to through-going muon ra- 



tio. Fewer upward stopping muons were observed than pre- 
dicted, while the observed number of upward through-going 
muons was consistent with the theoretical prediction within 
the errors. The observed ratio of stopping to through-going 
muons was 0.229 ± 0.015(ifar) ± 0.003(i3'.s), while the ex- 
pected ratio was 0.430 ±0.065. The expected ratio has the- 
oretical uncertainties from cross sections (±4.7 %), the cos- 
mic ray spectral index (±12.5 %), and the flux model depen- 
dence (±7.1 %). Fig. 1371 shows the expected ratio of stop- 
ping to through-going upward muon events as a function of 
Ain^ along with the measured ratio, which was smaller than 
the prediction by more than 3 standard deviations. The ob- 
served value can be explained assuming neutrino oscillations 
with Am^ in the range of 10^^ to 10^^ eV^. This stopping to 
through-going ratio is no longer explicitly used in the oscilla- 
tion fits, but is presented for comparison to older work 
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FIG. 34: Expected U /D ratio for FC single-ring jU-like + PC events 
as a function of Am^ for full ^ mixing. The data sample is 
divided into FC sub-GeV with < 400 MeV/c, FC sub-GeV with 
Pf, > 400 MeV/c, and FC multi-GeV + PC events. The ratio for the 
data together with the ±lcT statistical error are shown by the hor- 
izontal lines. The systematic errors are shown by the band in the 
expectation. 



VI. OSCILLATION ANALYSIS 

The observed deficits of muon neutrino interactions are in 
strong disagreement with the expectation in the absence of 
neutrino oscillations. Oscillation between electron neutrinos 
and muon neutrinos cannot explain the data, as no surplus of 
upward-going electron neutrinos is observed in the multi-Ge V 
data sample; an attempt at a two-flavor v^, ^ Vc fit results in 
a generally poor fit, with difference of more than 100 with 
respect to the ^ Vx analysis described below. A variety 
of exotic alternatives such as neutrino decay were considered, 
however, none fit the data as well as the ^ Vx scenario an- 
alyzed below. Atmospheric oscillation into Vx is mostly 
characterized by disappearance, as the majority of the flux 
is below the 3.5 GeV neutrino energy threshold for charged 
current i production. We carefully studied the alternative that 

could oscillate to a sterile neutrino state 1.8 1.1 . which would 
also result in disappearance. However, the lack of matter- 
induced suppression of oscillation and the relative up-down 
symmetry of the multi-ring sample with considerable neutral 
current fraction eliminated this hypothesis from serious con- 
siderations. The final Super-Kamiokande statistical analysis 
of these alternative scenarios, as well as the standard three 
flavor oscillation analysis, will appear in other publications. 
In this paper, we therefore establish the best-fit parameters of 
V/j ^ Vx oscillation. 

The analysis is based on a comparison between data and 
Monte Carlo, suitably binned to convey information about 
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FIG. 35: Event rates as a function of elapsed SK running time for 
upward-going muon data. 
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FIG. 36: Expected ratio of vertical to horizontal upward through- 
going muons as a function of Am^ for full ^ Vx mixing. Verti- 
cal and horizontal are defined to be the number of events in — 1 < 
cos < —0.5 and —0.5 < cos < 0.0, respectively. The ratio for the 
data together with ±la combined statistical and systematic error is 
also shown by the horizontal lines. The systematic error is shown by 
the band in the expectation. 



neutrino type, neutrino energy, and flight distance. The neu- 
trino type, Vf. or v^ is classified by the identification of the 
main Cherenkov pattern as showering or non-showering re- 
spectively. Penetrating particles such as upward-going muons 
and partially-contained events are assumed to arise from v^ 
interactions. The neutrino energy is correlated with the outgo- 
ing lepton momentum using the interaction models described 
in Section IV. The flight distance is correlated with the zenith 
angle as described by Figs. and |32] To study neutrino os- 
cillation using Eq.^ we reweight each simulated event using 
the Monte Carlo "truth" information of and L and bin the 
reweighted events for comparison with the detected data. Un- 
like our analysis using the ratio L/Zi 1 86], we make no attempt 
to estimate L or on an event-by-event basis." 
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FIG. 37: Expected ratio of stopping to tiirougli-going upward muons 
as a function of Am^ for full ^ mixing. Tiie ratio for the data 
together with ±la combined statistical and systematic error is also 
shown by the horizontal lines. The systematic error is shown by the 
band in the expectation. 

We used all of the data samples with a well-identified CC 
component, namely: FC single-ring /j-like, PC, multi-ring 
/j-like, upward stopping muons, and upward through-going 
muons. Because the flux of electron neutrinos provides a 
powerful constraint through the accurately predicted V^/Vg 
ratio, the single-ring e-like events were included in the fit. 
The FC single ring jU-like and e-like samples were divided in 
logarithmically-spaced momentum bins. All samples were di- 
vided in 10 zenith angle bins. In total 180 bins were used in 
the analysis: 150 for the FC sample, 10 for the PC sample, 10 
for the upward stopping muon sample, and 10 for the upward 
through-going muon sample. The number of observed and ex- 
pected events for each bin are summarized in the Appendix. 
A statistic is defined by the following sum: 

^" = L- -2 

A^f''P=A^f.p(v„^Vp) . (5) 

In the first sum, Nf^^ is the number of observed events in 
the /* bin and A^^^"'' is the expected number of events based 
on a Monte Carlo simulation and o, combines the statistical 
uncertainties in the data and Monte Carlo simulation. Dur- 
ing the fit, the values of A^J^"'' are recalculated to account for 
neutrino oscillations and systematic variations in the predicted 
rates due to uncertainties in the neutrino flux model, neutrino 
cross-section model, and detector response. is the number 
of events predicted from the MC without neutrino oscillation 
for the bin. The appearance of Vx as a result of oscillations 
is taken into account by adding into the Monte Carlo distri- 
butions simulated Vx interactions which pass all cuts. These 
events show up mainly in the multi-GeV e-like sample, but are 
not easily distinguished on an event-by-event basis. We are 
undertaking a separate analysis, to be published later, which 
will study Vx appearance in the atmospheric neutrino flux. 

The systematic uncertainties are represented by 39 param- 
eters 8;. During the fit, these 39 8; are varied to minimize 



for each choice of oscillation parameters sin^20 and Arn^. 
Among these, only 38 contribute to the x^, since the absolute 
normalization is allowed to be free. The factor /j represents 

the fractional change in the predicted event rate in the bin 
due to a variation of the parameter 8^. The second sum in the 

definition collects the contributions from the systematic un- 
certainties in the expected neutrino rates. The 8; are listed in 
Tables I VIII I Villi [Txl and IXI with their estimated uncertainties 
and the resulting best-fit values. Entries of the same number 
are treated as fully coiTelated although the effect of the uncer- 
tainty varies in size depending on its relative importance to the 
energy bin of certain sub-samples. For example, the source 
of the up/down uncertainty (No. 8) is due to the uncertainty 
in the geomagnetic field effect, especially above the Super- 
Kamiokande detector. The uncertainty is large for low energy 
neutrinos coming from primary cosmic rays below the geo- 
magnetic cutoff, but the effect of the uncertainty is decreased 
due to the large scattering angle in the neutrino interactions. 
As a result, events in the middle energy range are the most 
influenced by this particular systematic uncertainty. Refer to 
the footnotes in the tables for more detail. 

A global scan was made on a (sin^29,IogAm^) grid min- 
imizing at each point with respect to 39 parameters listed 
in Tables IVTII IVTni IPO and IXI At each grid point, the local 
minimum of x^ are derived by assuming a linear dependence 
of N^^'^ on each of the parameters. At the minimum %^ loca- 
tion, dx^/dEj = for each of the parameters 8/. As a result, 
the minimization of x^ in Eqn.l^is equivalent to solving the 
following system of k = 1 , 39 linear equations is^ll : 

39 r 1 180 / iuexP . f^^V , fi , fi\' 

180 ,')vobs_A,exp, exp w 

= Z - ' 2 ' ' (6) 

where Oj is the estimated uncertainty in the parameter e^-. One 
of Oj corresponds to the absolute normalization uncertainty. 
In this case, 1 /oj is set to 0, since the absolute normalization 
is a free parameter in our analysis. 

The minimum x^ value, xl,i„ = 174.8/177 DOF, is located 
at (sin2 2e = 1.00, Am^ = 2.1 x 10"^ eV^). The number 
of DOF is found by 180 terms in the %^ sum plus 38 sys- 
tamtic constraints in the x^ sum minus 39 minimized pa- 
rameters minus the two physics parameters of sin^29 and 
Ain^. The overall normalization is not used as a constraint 
to x^- The best-fit values of the parameters 8/ obtained at 
the global minimum are summarized in Tables I VIII IVIIII . Hxl 
and|X] For the most part, the parameters 8^ are fit within 
their estimated 1 a errors. Including the unphysical region 
(sin^29 > 1) in the scan, the minimum value is obtained 
at (sin2 2e 1.02,Am2 = 2.1 x 10"^ eV^). The minimum 
X^ value, Xmin ~ 174.5/177 DOF, in the unphysical region is 
lower than that in the physical region by 0.29 . Contours cor- 
responding to the 68%, 90% and 99% confidence intervals are 
located at xL«+ 2.60, 4.98, and 9.60 respectively, where Xmi„ 
is the minimum x^ value in the physical region and are shown 
in Fig.|38l These intervals are derived based on a two dimen- 
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a(%) best-fit No. 

(A) Systematic uncertainties in neutrino flux 



Absolute normalization 


free 


11.9 


1 




£v < 5 GeV 


3.0 


-2.4 


2 




£-v > 5 GeV 


3.0* 


0.1 


3 




£v < 10 GeV 


5.0 


1.5 


4 




£v > 10 GeV 


5.0^^ 


0.0 


5 




£v < 10 GeV 


5.0 


-1.3 


6 




£v > 10 GeV 


5.0*" 


0.9 


7 


Up/dowr/ 


<400MeV e-like 


0.5 


0.2 


8 




jU-like 


0.8 


0.3 


8 




>400MeV e-like 


2.1 


0.9 


8 




ii/-like 


1.8 


0.8 


8 




Multi-GeV e-like 


1.5 


0.7 


8 




jU-like 


0.8 


0.3 


8 




PC 


0.4 


0.2 


8 




Sub-GeV multi-ring fA 


0.8 


0.3 


8 




Multi-GeV multi-ring fj 


0.7 


0.3 


8 


Horizontal/vertica/ 


<400MeV e-like 


0.3 


0.0 


9 




jU-like 


0.3 


0.0 


9 




>400MeV e-like 


1.2 


0.1 


9 




iU-like 


1.2 


0.1 


9 




Multi-GeV e-like 


2.8 


0.2 


9 




/i-like 


1.9 


0.1 


9 




PC 


1.4 


0.1 


9 




Sub-GeV multi-ring fj 


1.5 


0.1 


9 




Multi-GeV multi-ring jj 


1.3 


0.1 


9 


K/tz ratio « 




20.0 


-6.3 


10 


Ly (production height) 


lO.O'' 


-0.6 


11 


Energy spectrum' 


Ek < 100 GeV 


0.03 


0.031 


12 




Ek > 100 GeV 


0.05 


0.052 


12 


S ample-by- sample' 


FC Multi-GeV 


5.0 


-5.2 


13 




PC + upward stopping /j 


5.0 


-3.9 


14 



"A positive number means the number of MC + V,, events is increased. 

''Error linearly increases with log£v from 3 % (5 GeV) to 10 %(100GeV). 

''A positive number means the number of MC (v^) events is increased. 

''Error linearly increases with log£v from 5 %(10 GeV) to 10 %(100 GeV). 

""Error lineai'ly increases with log£v from 5 %(10 GeV) to 25 %(100 GeV). 

^Up/d own (horizontal/vertical) uncertainty in neutrino flux is assumed to be 
fully correlated. All of the samples listed are simultaneously vaiied according 
to the systematic uncertainty factors. A positive number means the number 
of MC upward (horizontally-going) events is increased. 

*20 % uncertainty in K/n production ratio in cosmic ray interactions in the 
atmosphere. A positive number means that the fraction of K is increased. 

''10% uncertainty in the atmospheric density stiucture. A positive number 
means a more compressed atmospheric density structure. 

'0.03 and 0.05 uncertainties in the spectral index of the primary cosmic 
rays below and above 100 GeV. Spectral index uncertainties below and above 
100 GeV are assumed to be correlated. A positive number means that the 
spectrum is harder The predicted flux was changed ai'ound an arbitrary ref- 
erence energy of 10 GeV. 

^Different flux calculations predict different energy dependences that cannot 
be explained by a simple spectral index uncertainty. See the lower panel 
of Fig.|2] From a comparison of the predicted number of events based on 
different flux models, 5 % is assigned as the relative normalization uncertainty 
for these samples. 

TABLE VII: Summary of systematic uncertainties in the prediction 
of the atmospheric neutrino flux. Estimated uncertainty and the best- 
fit value are listed for each error. The last column shows the error 
parameter numbers (_/), which appeared in Eqs|4]and|6| 



o(%) best-fit No. 



(B) Systematic uncertainties in neutrino interaction 



Ma in quasi-elastic and single-n 


10.0" 


0.5 


15 


Quasi-elastic scattering (model dependence) 


I.O* 


-0.95 


16 


Quasi-elastic scattering (cross-section) 


lO.O 


5.6 


17 


Single-meson production (cross-section) 


lO.O 


-4.7 


18 


A'lulti-pion production (model dependence) 


I.O*^ 


1.47 


19 


Multi-pion production (total cross-section) 


5.0 


-0.2 


20 


Coherent pion production (total cross-section) 


30.0 


0.4 


21 


NC/CC ratio 


20.0 


2.9 


22 


Nuclear effect in "'0 '^ 


30.0 


-7.2 


23 


Energy spectrum of pions 


1.0^ 


0.50 


24 


CC Vt interaction cross section 


30.0 


0.2 


25 



"10 % uncertainty in the axial vector mass. Ma (See Sec. lIVi . value. 

''Difference from the model in Ref. |82ll is set to 1.0 . 

"^Difference from the model in Ref. I83ll is set to 1.0 . 

''a positive number means more NC events in the Monte Carlo. 

""30 % uncertainty in the mean free path of hadrons in the '^O nucleus. A 
positive number means stronger nuclear effect in '^O. 

'^The difference in the predicted pion energy spectrum by NEUT and 
NUANCE interaction models is taken as 1 standard deviation, and is set to 
1.0. 

TABLE VIII: Summary of systematic uncertainties in neutrino in- 
teractions. Estimated uncertainty and the best-fit value are listed for 
each error. The last column shows the error parameter numbers (j), 
which appeared in Eqs|4|and|6| 



sional extension of the method described in Ref. jssll . Fig- 
ure |39l shows the —Xmin distributions projected to sin^29 
and Am^ axes, in which the minimum — Xmin values for 
each sin^ 29 and Am^ are plotted. The x^ ~ xLn distribution 
is rather flat between Am^ = 2. Ox lO^^eV^ and 2.5 x lO^^eV^. 
Any Am^ in this range fits the data nearly as well as the best-fit 
point. 

Assuming no oscillation, (sin^ 29 = 0, Am^ = 0), we found 
a x^ value of 478.7 for 179 DOF, where only the overall nor- 
malization is a free parameter. We allowed all systematic un- 
certainty terms to be minimized, yet the fit was greatly inferior 
to the best-fit including neutrino oscillations. 

We have also estimated the allowed neutrino oscillation pa- 
rameters by performing the same fitting procedure using in- 
dependent subsamples of the data: FC single-ring sub-GeV 
events below 400 MeV/c, FC single-ring sub-GeV events 
above 400 MeV/c, FC single-ring multi-GeV events, PC 
events, FC multi-ring events, and upward-going muon events. 
In each independent fit, only the relevant parameters out of the 
set of 39 were minimized. The results are shown in Fig. 1401 
The allowed region contours found by fitting these six sub- 
samples are consistent with each other and with the combined 
fit to all events. 

In addition, the same oscillation analyses were repeated us- 
ing different flux models (but with the same neutrino interac- 
tion Monte Carlo program) and different neutrino interaction 
Monte Carlo program (but with the same flux model). The 
90 % C.L. allowed parameter regions are compared in Fig. 1411 
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a (%) best-fit No. 
(C) Systematic uncertainties in event selection 



Reduction for fully-contained event 


0.2 


0.0 


26 


Reduction for partially-contained event 


2.6 


0.3 


27 


Detection efficiency" 


upward stopping /j 


1.3 


-0.2 


28 




upward through-going /j 


0.5 


-0.1 


28 


FC/PC separation* 




0.9 


-0.3 


29 


Hadron simulation 




1.0'' 


-0.24 


30 


Non-v BG'' 


Sub-GeV e-like 


0.4 


0.1 


31 




/j-like 


0.1 


0.0 


32 




Multi-GeV e-like 


0.2 


0.0 


31 




fj-like 


0.1 


0.0 


32 




PC 


0.2 


0.0 


32 


Upward stopping/through-going /j separation " 


0.4 


0.0 


33 



"Goodness of upward-going fit is used to select the upward-going /j sam- 
ple. The difference of the goodness between the data and MC is considered 
as the source of the uncertainty in the detection efficiency. Uncertainties for 
upward stopping /j and upward through-going are assumed to be correlated. 

''The number of hits in the OD cluster is used to separate the FC and PC 
events. See Fig. 1161 The systematic uncertainty in the number of hits in the 
OD cluster causes 0.9 % uncertainty in the number of the PC events. The 
number of FC events changes anti-correlated with the change in the number 
of PC events. A positive number means that the number of MC FC events is 
increased. 

''Difference from the FLUKA model. A positive number means more 
hadrons, mostly pions, in neutral current interactions are identified as ;U-like. 

''The background sources are flasher PMTs and neutron interactions for e- 
hke events and cosmic ray muons for ^i-like events. It is assumed that the 
background sources are un-correlated between e-like and /.i-like events. The 
background for sub- and multi-Ge'V samples in the e-like and ;U-like events 
ai'e assumed to be correlated. The background for the PC sample is also 
assumed to be correlated with the FC ;U-like samples. Only positive numbers 
are allowed for the background. 

''The number of hits in the OD cluster at the exit point of a muon is used 
to separate the upward stopping and through-going muon events. The uncer- 
tainty in the number of hits in the OD cluster causes 0.4 % uncertainty in the 
stopping/through-going ratio. A positive number means that the number of 
MC stopping muons is increased. 

TABLE IX: Summary of systematic uncertainties in event selection. 
Estimated uncertainty and the best-fit value are listed for each er- 
ror. The last column shows the error parameter numbers (j), which 
appeared in Eqs|4|and|6| 



O (%) best-fit No. 
(D) Systematic uncertainties in event reconstruction 



JVlIlg aCUtli d.Li(JIl \. H-UL/ iVlC V 




U.J 


9 6 


34 




/j-like 


2.4 




34 


>400MeV 


e-like 


3.4 


1.4 


34 




/j-like 




n 5 
vj.j 


34 


MuIti-GeV 


e-like 




U.J 


34 




/j-like 


6.2 


2.5 


34 


Sub-GeV multi-rinj 




^ 7 
J. 1 


-1 S 


34 


Multi-GeV multi-ring jj 


7.2 


-7 


34 


Particle identification* Sub-GeV 


e-like 


n 6 


u.z, 


35 




/J-like 


0.6 


-0.2 


35 


Multi-GeV 


e-like 


0.4 


0.1 


35 




jj-Yike. 


0.4 


-0.1 


35 


Sub-GeV multi-rinj 




3.4 


-0.9 


36 


Multi-GeV multi-ring jj 


4.7 


-1.2 


36 


Energy calibration for FC event ' 




2.0 


0.4 


37 


Energy cut for upward stopping muon 




1.1 


-0.2 


38 


Up/down symmetry of energy calibration ^ 




0.6 


0.0 


39 



"Ring separation uncertainty is assumed to be fully correlated. Namely, 
if the number of single-ring sub-GeV e-like events have to be increased, the 
number of single-ring multi-GeV e-like events and single-ring sub- and multi- 
GeV fj-liks events have to be increased according to the systematic uncer- 
tainty factors. On the other hand, in this case, the number of multi-ring /(-like 
events have to be decreased. A positive number means the number of MC 
events for the corresponding sample is increased. 

''The particle identification uncertainty is anti-correlated between e-like and 
/.i-like events. It is assumed that the particle identification uncertainty is cor- 
related between sub- and multi-GeV energy regions. However, it is assumed 
that it is not correlated between single- and multi-ring events. A positive 
number means the number of MC events for the corresponding sample is in- 
creased. 

'2 % uncertainty in the absolute energy scale of the detector A positive 
number coiTesponds to increasing the visible energy of MC events. 

''a positive number means that the energy of MC events is increased for 
upward-going direction. 

TABLE X: Summary of systematic uncertainties in event reconstruc- 
tion. Estimated uncertainty and the best-fit value are listed for each 
error. The last column shows the error parameter numbers (j), which 
appeared in Eqs|4|and|6| 



The allowed regions from these analyses overlap well, demon- 
strating that the measured parameters do not strongly depend 
on the choice of flux or interaction model from which we 
start the fitting procedure. However, the allowed region ob- 
tained based on the flux model of Ref. allows for slightly 
higher Arn^. We studied the reason for this difference in detail, 
and found that the main reason was the slightly harder energy 
spectrum in the upward-going muon energy range (Fig.|2ll. 

Finally, we point out that a separate L/E analysis of the 
same running period |86], using only selected high resolution 
FC and PC events, gave an allowed oscillation parameter re- 
gion consistent with this result. This is shown in Fig.|^ with 
a magnified view of the region and a linear scale in Am^. The 
L/E analysis provided a slightly better constraint in Arn^ due 



to locating the oscillatory dip; the present analysis constrains 
sin^ 29 better due to high statistics in the up-down asymmetry. 



VII. CONCLUSION 

Super-Kamiokande has observed more than 15,000 atmo- 
spheric neutrino events during the first data taking period 
between 1996 and 2001. Atmospheric neutrino events ob- 
served in Super-Kamiokande have an energy range from 
about 100 MeV to 10 TeV, and a neutrino flight-length from 
about 10 km to 13,000 km. These wide energy and flight- 
length ranges together with high statistics made it possible to 
study neutrino oscillations. Especially, the predicted up-down 
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FIG. 38: Allowed oscillation parameters for ^ Vx oscillations. 
Three contours correspond to the 68% (dotted line), 90% (solid line) 
and 99% (dashed line) C.L. allowed regions. 



FIG. 40: 90 % confidence level allowed oscillation parameter regions 
for ^ Vt oscillations from six sub-samples. In this plot, 90 % 
confidence interval is defined to be = Xm/n +4.61, where x^„„ is 
the minimum jj- value including the unphysical parameter region. 
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FIG. 41: Left: 90% confidence level allowed oscillation parameter 
regions for ^ Vt oscillations, based on the NEUT neutrino interac- 
tion model, from different flux models (solid line; |2^ . dashed line; 
E], dotted line; |29]). Right: The 90% C.L. allowed regions based 
on a different neutrino interaction model (NUANCE |48]) for FC-l-PC 
events with the flux model of Ref. 1 28] (dashed line) is compared 
with that based on NEUT with the same flux. In this plot, Monte 
Carlo events from NEUT were used for upward-going muons. 



FIG. 39: ^Xmin projected onto the sin^2e and Ain^ axes. The 
minimum value at each sin^ 26 and Arn^ is plotted. 



asymmetry of the atmospheric neutrino flux enabled us to ac- 
curately estimate the mixing parameter sin^ 29. The observed 
muon neutrino events showed a clear zenith angle and energy 
dependent deficit of events, while the electron neutrino events 



were in good agreement with the prediction. A detailed neu- 
trino oscillation analysis confirmed that the full data set was 
explained well by ^ Vx oscillations. Various systematic 
effects were included in the oscillation analysis. The mea- 
sured neutrino oscillation parameters were sin^ 29 > 0.92 and 
1.5 X 10-3 < An? < 3.4 x lO^^ eV^ at 90% C.L. This result 
gives the most accurate determination of sin^ 29 and is consis- 
tent with the somewhat more accurate measurement of Arn^ 
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FIG. 42: The 68, 90 and 99 % confidence level allowed oscilla- 
tion parameter regions obtained by an L/E analysis and by the 
present analysis are compared. 



determined by the independent study of high resolution L/E 
events I861 . 
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A. Appendix 

Table lXII summarizes the number of observed and expected 
FC and PC events for each bin. The Monte Carlo prediction 
does not include neutrino oscillations. Table IXIll summarizes 
those for upward-going muons. These binned data are used in 
the oscillation analysis. Table IXIIll summarizes neutrino en- 
ergy at which 50% of events are accumulated for each energy 
bin in the absence of neutrino oscillations. The fraction of 
various neutrino interaction modes are also listed. 
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FC single-ring e-like 





I 


II 


m 


IV 


V 


VI 


vn 


nx 


IX 


X 


1 

2 
3 
4 
5 
6 
7 


114(79.29) 
96(75.62) 
76(64.16) 
48(45.35) 
26(21.68) 
33(29.29) 
23(21.36) 


95(83.33) 
93(71.70) 
80(66.93) 
57(47.92) 
35(23.19) 
35(33.22) 
31(25.66) 


74(81.41) 
96(73.22) 
80(65.78) 
62(50.12) 
31(25.13) 
41(34.92) 
28(30.10) 


94(82.04) 
90(69.44) 
69(63.60) 
52(50.91) 
37(25.84) 
37(39.72) 
42(37.93) 


88(83.99) 
89(68.36) 
72(64.57) 
60(51.62) 
24(25.55) 
46(42.84) 
63(45.70) 


91(79.78) 
85(68.75) 
60(64.06) 
74(51.60) 
38(25.93) 
49(43.86) 
37(42.99) 


79(79.49) 
85(69.49) 
69(62.39) 
55(50.75) 
34(24.96) 
49(40.65) 
54(35.42) 


74(84.19) 
74(67.18) 
71(61.66) 
58(49.13) 
24(26.14) 
32(39.50) 
34(32.54) 


91(81.50) 
83(71.14) 
85(59.72) 
60(46.51) 
21(23.63) 
36(32.03) 
22(26.00) 


100(82.93) 
78(69.68) 
63(57.48) 
43(42.45) 
20(18.51) 
36(27.30) 
18(19.48) 


FC single-ring /7-like 


1 

2 
3 
4 
5 

6+7 


36(54.73) 
86(124.32) 
94(118.74) 
52(91.07) 
27(43.35) 
34(89.46) 


40(53.66) 
77(122.76) 
60(112.36) 
48(87.99) 
22(45.89) 
46(86.45) 


39(54.39) 
99(123.28) 
81(113.06) 
53(90.52) 
22(44.91) 
42(86.48) 


37(55.09) 
86(121.54) 
94(115.77) 
53(91.04) 
37(44.48) 
49(92.54) 


35(55.75) 
87(119.12) 
87(112.80) 
68(94.68) 
25(47.03) 
54(96.41) 


34(53.79) 
80(119.67) 
84(112.77) 
68(91.11) 
40(47.52) 
73(94.68) 


35(53.46) 
91(122.64) 
116(113.40) 
72(89.57) 
41(47.91) 
95(96.23) 


45(53.57) 
85(117.84) 
119(111.91) 
81(88.15) 
41(42.59) 
87(88.39) 


48(52.63) 
94(115.67) 
97(108.09) 
91(84.45) 
46(44.04) 
78(84.46) 


46(52.10) 
76(120.59) 
118(104.54) 
86(82.90) 
48(43.60) 
93(84.82) 


FC multi-ring /i-like 


sub-GeV 
multi-GeV 


14(27.57) 
18(63.27) 


8(31.15) 
29(63.34) 


20(33.42) 
31(66.89) 


14(33.68) 
28(75.69) 


25(36.13) 
41(86.27) 


16(35.59) 
69(82.70) 


21(34.12) 
55(77.10) 


32(32.94) 
54(73.27) 


29(28.85) 
59(62.44) 


29(29.12) 
55(60.87) 


PC 




49(88.97) 


45(88.81) 


59(104.46) 


89(125.12) 


117(154.61) 


156(158.34) 


114(128.39) 


109(103.00) 


85(89.70) 


88(88.18) 



TABLE XI: Summary of the number of observed (MC expected) FC and PC events for each bin. Neutrino oscillation is not included in the 
Monte Carlo prediction. Roman numbers represent zenith angle regions equally spaced between cos® = — 1 and cos® = 1. The numbers 

in the log{Pigp) column show the momentum ranges. The momentum ranges are <250, 250-400, 400-630, 630-1000 and >1000 MeV/c in 
sub-GeV samples for momentum range numbers 1 to 5 and <2500 and >2500 MeV/c in multi-GeV samples for momentum range numbers 6 
to 7. 



upward through-going muon 





I 


II 


III 


IV 


V 


VI 


VII 


IIX 


IX 


X 


# Events 


85 


113 


116 


138 


159 


183 


178 


267 


286 


316.6 


# Expected 


96.13 


114.85 


122.48 


136.88 


145.79 


169.17 


187.32 


210.92 


228.76 


257.16 


Efficiency 


95.2% 


93.9% 


92.4% 


95.9% 


94.0% 


97.2% 


96.9% 


99.0% 


96.2% 


95.5% 


Flux 


0.862 


1.060 


1.045 


1.216 


1.388 


1.589 


1.557 


2.365 


2.574 


2.953 


Stat. err. 


±0.093 


±0.100 


±0.097 


±0.104 


±0.110 


±0.117 


±0.117 


±0.145 


±0.152 


±0.191 


Expected 


0.976 


1.078 


1.103 


1.206 


1.274 


1.469 


1.638 


1.868 


2.060 


2.393 


Stat. err. 


±0.021 


±0.021 


±0.021 


±0.022 


±0.022 


±0.024 


±0.025 


±0.027 


±0.029 


±0.032 










upward stopping 


muon 










# Events 


28 


23 


37 


30 


27 


37 


37 


48 


65 


85.7 


# Expected 


51.24 


54.06 


56.69 


64.96 


67.60 


68.21 


78.85 


80.96 


94.00 


96.90 


Efficiency 


99.8% 


99.1% 


99.8% 


108.1% 


102.1% 


101.2% 


103.5% 


103.1% 


105.0% 


100.2% 


Flux 


0.286 


0.217 


0.337 


0.265 


0.236 


0.322 


0.323 


0.425 


0.589 


0.807 


Stat, err 


±0.054 


±0.045 


±0.055 


±0.048 


±0.045 


±0.053 


±0.053 


±0.061 


±0.073 


±0.172 


Expected 


0.523 


0.509 


0.517 


0.574 


0.591 


0.594 


0.689 


0.715 


0.851 


0.913 


Stat. err. 


±0.015 


±0.015 


±0.015 


±0.015 


±0.015 


±0.015 


±0.016 


±0.017 


±0.019 


±0.020 



TABLE XII: Summary of the number of observed and expected upward-going muons for each bin, efficiencies, and the corresponding 
flux. Neutrino oscillation is not included in the Monte Carlo prediction. The errors on the observed fluxes are statistical, the units of flux 
X 10~'^cm~'^s~^sr~^. The Roman numerals refer to zenith angle regions equally spaced between cos® = —1 and cos® = 0. 
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FC single-ring e-like 





fvCGeV) 


CCv,(%) 


CC (%) 


NC (%) 


1 


sub-GeV 100-250 


0.31 


87.6 


2.2 


10.1 


2 


250-400 


0.48 


89.1 


1.5 


9.5 


3 


400-630 


0.72 


88.7 


1.8 


9.5 


4 


630-1000 


1.1 


86.8 


3.2 


10.1 


5 


>1000 


1.5 


86.5 


4.5 


8.9 


6 


multi-GeV <2500 


2.3 


85.7 


5.8 


8.5 


7 


>2500 


5.4 


79.0 


8.4 


12.6 




FC sinj 


jle-ring jU-like 






1 


sub-GeV 200-250 


0.50 


0.3 


90.9 


8.8 


2 


250-400 


0.68 


1.2 


95.6 


3.3 


3 


400-630 


0.86 


0.7 


97.7 


1.6 


4 


630-1000 


1.2 


0.5 


99.0 


0.4 


5 


>1000 


1.5 


0.4 


99.3 


0.2 


6+7 


multi-GeV 


2.6 


0.4 


99.4 


0.2 


FC multi-ring /j-like 




sub-GeV 


1.9 


3.6 


90.5 


5.9 




multi-GeV 


3.6 


2.3 


94.9 


2.7 


PC 




7.9 


1.8 


97.3 


0.9 


Upwarc 


-going muons 




upward stopping muon 


11.1 


1.0 


98.6 


0.4 




upward through-going muon 


113.5 


0.2 


99.7 


0.1 



TABLE XIII: Description of momentum bins used for this analysis, corresponding to raws of Tables IXTI and IXlTI Also tabulated are the 
medium parent neutrino energy, and the relative fractions of CC Ve, CC and NC interactions in the absence of neutrino oscillations as 
estimated by the Monte Carlo program. 
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